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NOMENCLATURE 


R<* 

=  real  vector  space  of  dimension  a  ^ 
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=  belongs  to 

A  [  *  ] 

=  Eigenvalues  of  the  matrix  [ - ] 

cr  [  -  ] 

=  singular  value  of  the  matrix  [•] 

=  {A([-][']t)}V2 
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=  symmetric  part  of  a  matrix  [ • ] 

K-)| 

=  modulus  of  the  entry  ( • ) 

[•In 

=  modulus  matrix  =  matrix  with  modulus  entries 

V 

=  for  all 

xa 

=  a  x  a  Identity  matrix 

IJC 

=  International  Journal  of  Control 

IEEE  AC 

=  IEEE  Transactions  on  Automatic  Control 
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INTRODUCTION  AND  PERSPECTIVE 


I. 


It  is  well  known  that  the  inaccuracies  in  the  mathematical  models  of 
physical  systems  can  severely  compromise  the  resulting  control  designs.  The 
errors  associated  with  mathematical  models  of  physical  systems  may  be  broadly 
categorized  as  i)  parameter  errors,  ii)  truncated  models  (errors  in  model 
order),  iii)  neglected  or  incorrectly  modeled  external  distrubances  and  iv) 
neglected  nonlinearities.  It  is  the  inevitable  presence  of  these  errors  in 
the  model  used  for  design  that  eventually  limits  the  performance  attainable 
from  the  control  system  designs  produced  by  either  classical  (frequency  domain) 
or  modern  (time  domain)  control  theory.  The  problem  of  model  errors  is  more 
critical ,  in  general,  for  large  scale  Linear  Quadratic  Gaussian  (LQG)  optimal 
control  problems  and  in  particular,  for  aerospace  applications  like  Large 
Space  Structure  (LSS)  control  and  other  aeroelastic  systems.  These 
applications  are,  of  course,  of  extreme  importance  to  the  U.S.  Air  Force.  One 
fundamental  problem  of  LSS  control  is  the  control  of  a  large  dimensional 
system  with  a  controller  of  much  smaller  dimension  (model /controller 
truncation)  compounded  with  modal  data  uncertainty  (parameter  errors).  In  the 
light  of  these  observations,  it  is  evident  that  Robustness1  is  an  extremely 
desirable  (sometimes,  necessary)  feature  of  any  feedback  control  design 
proposed  for  LSS  control.  'Robustness'  studies  of  LQG  regulators  is  the 
central  theme,  of  the  present  research. 

For  our  present  purposes  a  'robust'  control  design  is  that  design  which 
behaves  in  an  'acceptable'  fashion  (i.e.  satisfactorily  meets  the  system 
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specifications)  even  in  the  presence  of  modeling  errors.  Since  the  system 
specifications  could  be  in  terms  of  stability  and/or  performance  (regulation, 
time  response,  etc.)  we  can  conceive  two  types  of  robustness,  namely, 
'Stability  Robustness'  and  'Performance  Robustness*.  Limiting  our  attention 
in  this  research  to  'parameter  errors'  and  'mode  truncation'  as  the  two  types 
of  modeling  errors  that  may  cause  instability  (or  performance  degradation)  in 
the  system,  we  formally  define  'stability  robustness'  and  'performance 
robustness'  as  follows: 

'Stability  Robustness1;  Maintaining  closed  loop  system  stability  in  the 
presence  of  modeling  errors  mainly  parameter  variations  and  model/controller 
truncation . 

'Performance  Robustness':  Maintaining  satisfactory  level  of  performance  (or 
regulation)  in  the  presence  of  modeling  errors  mainly  parameter  variations  and 
mode  truncation. 

Clearly,  the  definition  of  'Performance  (or  Regulation)  Robustness* 
embodies  the  definition  of  'Stability  Robustness'.  If  the  main  concern  is  for 
stability  alone,  then  one  can  analyze  and  synthetize  for  robust  stability 
whereas  if  both  stability  and  performance  (regulation,  for  the  case  of  linear 
regulators)  are  of  concern,  then  one  needs  to  resort  to  the  'Performance 
Robustness'  aspect.  In  this  report  these  two  aspects  are  clearly  delineated 
and  accordingly  are  given  separate  treatment. 

The  published  literature  on  the  'robustness*  of  linear  systems  can  be 
viewed  mainly  from  two  perspectives,  namely  i)  frequency  domain  analysis  and 
ii)  time  domain  analysis.  The  main  direction  of  research  in  frequency  domain 
has  been  to  extend  and  generalize  the  well  known  classical  single  input 
single  output  treatment  to  the  case  of  multiple  input  multiple  output  systems, 
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using  the  singular  value  decomposition.  In  the  case  of  frequency  domain 
results,  the  perturbations  are  mainly  viewed  in  terms  of  'gain1  and  'phase* 
changes.  The  time  domain  treatment  is  more  or  less  similar  to  the  frequency 
domain  treatment  in  spirit  but  quite  different  in  detail.  The  time  domain 
treatment  is  more  amenable  to  treating  perturbations  in  the  form  of  parameter 
variations,  nonlinearities  and  external  disturbances  and  also  for  the  physical 
interpretation  of  many  real  life  perturbations.  This  research  treats  the 
robustness  analysis  and  design  from  time  domain  viewpoint  and  in  particular 
focuses  on  the  well  known  Linear  Quadratic  Regulator  problem.  In  addition, 
the  main  tool  used  is  the  Lyapunov  stability  analysis  which  allows  time 
varying  perturbations  to  be  considered  in  the  analysis. 

The  problem  of  maintaining  the  stability  of  a  nominally  stable  system 
subject  to  perturbations  has  been  an  active  topic  of  research  for  quite  some 
time.  One  factor  which  clearly  influences  this  type  of  analysis  is  the 
characterization  or  type  of  'perturbation*.  Even  in  the  context  of  nominally 
stable  linear  systems,  the  'perturbations'  can  take  different  forms  like 
linear,  nonlinear,  time  invariant,  time  varying,  structured  and  unstructured. 
Structured  perturbations  are  those  for  which  bounds  on  the  individual  elements 
of  the  perturbation  matrix  are  known  (or  derived)  whereas  unstructured 
perturbations  are  those  for  which  only  a  norm  bound  on  the  perturbation  matrix 
is  known  (or  derived).  In  this  research,  we  focus  our  attention  on  linear, 
time  varying,  structured  perturbations  as  affecting  a  nominally  stable  linear 
time  invariant  system. 

With  this  perspective  in  mind,  the  report  is  organized  as  follows: 

Section  II  presents  the  analysis  for  stability  robustness  of  linear  systems. 
After  a  brief  literature  search,  upper  bounds  on  the  linear,  time  varying. 
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structured  perturbation  of  an  asymptotically  stable  linear  time  invariant 
system  are  presented  to  maintain  stability.  Then  a  state  transformation 
technique  is  presented  to  further  reduce  the  conservatism  of  the  Lyapunov 
stability  bound.  The  analysis  is  then  applied  to  ultimate  boundedness  control, 
the  analysis  of  which  also  depends  on  Lyapunov  stability  bound.  The 
usefulness  of  the  perturbation  bound  for  structured  uncertainty  is 
demonstrated  by  extending  the  analysis  to  present  sufficient  conditions  for 
the  stability  of  an  interval  matrix.  Section  III  is  completely  devoted  to  the 
design  of  controllers  for  robust  stability  along  with  an  aircraft  control 
example.  Then  in  Section  IV  attention  is  focused  on  the  performance 
robustness  analysis,  which  presents  upper  bounds  for  robust  regulation  and 
their  relationship  to  the  bounds  for  robust  stability.  Section  V  presents  a 
procedure  for  designing  controllers  for  robust  regulation  along  with  a  simple 
example  to  illustrate  the  procedure.  Section  VI  addresses  the  application  of 
this  'Perturbation  Bound  Analysis'  to  Large  Space  Structures  (LSS)  models  in 
which  uncertain  modal  data  and  mode  truncation  are  modeled  as  additive 
perturbations  and  the  results,  obtained  by  exploiting  the  special  structure  of 
LSS  models,  are  discussed.  Finally  Section  VII  offers  some  concluding  remarks 
and  explores  avenues  for  future  research  that  needs  the  continued  sponsorship 
of  the  Air  Force. 
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II.  ANALYSIS  FOR  STABILITY  ROBUSTNESS  OF  LINAER  SYSTEMS 


In  the  present  day  applications  of  linear  systems  theory  and  practice, 
one  of  the  challenges  the  designer  is  faced  with  is f  to  be  able  to  guarantee 
'acceptable'  behavior  of  the  system  even  in  the  presence  of  perturbations. 

The  fundamental  'acceptable'  behavior  of  any  control  design  for  linear  systems 
is  'stability'  and  accordingly  one  of  the  important  tasks  of  the  designer  is 
to  assure  'stability'  of  the  system  subject  to  perturbations. 

In  particular,  as  discussed  in  the  introduction,  we  concentrate  on 
'parameter  uncertainty'  as  the  type  of  perturbation  acting  on  the  system. 

This  section,  thus,  addresses  the  analysis  of  'stability  robustness'  of  linear 
systems  subject  to  parameter  uncertainty. 

2_^1_ _ Brief  Review  of  Literature 

The  aspect  of  robust  control  design  for  linear  regulators  subject  to 
parameter  uncertainty  has  also  been  an  active  topic  of  research  in  recent 
years.  Much  of  this  research  in  frequency  domain  is  carried  out  using 
singular  value  analysis.  The  special  issue  [1],  [2]  contain  valuable 

references  in  this  category.  A  recent  contribution  in  the  frequency  domain  is 
the  work  of  Daniels  and  Kouvaritakis  [3]  who  consider  'stability  margin'  both 
from  tolerance  to  uncertainty  as  well  as  typical  gain,  phase  margin  points  of 
view.  On  the  other  hand,  there  is  equally  interesting  research  reported  in 
time  domain.  Davison  [4]  focuses  on  the  structure  or  characterization  of  the 
robust  controller  and  presents  necessary  and  sufficient  conditions  for  such  an 
observer  based  controller  to  exist.  It  is  assumed  that  the  system  is  stable 
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under  the  given  range  of  parameter  perturbations  and  no  explicit  bounds  on  the 


tolerable  perturbations  for  stability  are  given.  The  parameter  space  design 
of  Ackermann  [5]  treats  only  single  input  systems,  with  time  invariant 
perturbations,  whereas  this  report  considers  time  varying  uncertain  parameters 
which  vary  within  certain  bounds.  Lunze  [6],  using  comparison  systems 
concept  and  Owens  and  Chotai  [7]  by  combining  frequency  domain  and  time 
domain  treatments,  present  attractive  schemes  involving  more  general 
perturbations.  Also  Desoer  et  al  [8]  have  established  conditions  for 
stability  robustness  of  linear  multivarible  interconnected  systems  for 
sufficiently  small  perturbations. 

From  a  different  perspective,  some  researchers  have  presented  analysis 
and  design  procedures  for  tolerable  perturbtions  for  robust  stability.  Kantor 
and  Andres  [9]  present  an  algorithm  to  determine  tolerable  perturbations  in 
frequency  domain  using  M  matrix  analysis.  In  time  domain  Horisberger  and 
Belanger  [10]  present  an  algorithm  to  determine  an  output  feedback  control 
gain  that  yields  the  largest  possible  tolerable  perturbation  such  that  the 
closed  loop  system  is  stable,  but  no  explicit  bounds  are  given.  Zheng  [11] 
presents  a  procedure  to  find  the  stability  regions  as  a  function  of  parameters 
but  considers  only  time  invariant  perturbations  and  no  synthesis  procedure  is 
given.  Eslami  and  Russell  [12]  also  address  the  same  problem  but  no  explicit 
bounds  are  obtained.  Chang  and  Peng  [13],  Patel  and  Toda  [14],  Patel,  Toda 
and  Sridhar  [15]  give  explicit  bounds  on  the  linear,  time  varying  perturbtions 
of  an  asymptotically  stable  linear  time  invariant  system  such  that  stability 
is  guaranteed  using  Lyapunov  theory.  However  the  same  measure  is  used  to 
propose  bounds  for  both  'structured*  as  well  as  ’unstructured*  perturbations 
[16],  making  the  bound  for  structured  perturbation  more  conservative. 
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In  what  follows  an  analysis,  carried  out  in  the  time  domain,  for  both 
types  of  perturbation  is  presented.  For  structured  perturbations,  an  improved 
upper  bound  is  presented  (over  that  of  [14]),  the  expression  for  which  is  such 
that  it  garners  the  full  'structural*  information  about  the  nominal  (and 
consequently,  the  perturbation)  matrix  into  it.  For  unstructured  perturbation 
case,  some  special  cases  of  the  nominally  stable  matrix  are  considered  for 
which  it  is  possible  to  give  simplified  upper  bounds  without  the  need  to  solve 
the  Lyapunov  equation.  Examples  presented  include  comparisons  with  the 
existing  aproaches. 

jy2* 2  Perturbation  Bounds L Jfor_J^inear  State  Space  Models 
Let  us  consider  the  following  linear  dynamic  system, 

x(t)  =  A(t)  x ( t )  (2.1a) 

=  [Aq  +  E(t>]  x(t)  (2.1b) 

A0  is  an  nxn  nominally  asymptotically  stable  matrix  and  E  is  an  nxn  'Error' 
matrix.  In  the  case  of  'structured'  perturbation,  the  elements  of  E  are  such 
that 

A  A 

Eij  (t)  <  eij  =  max  I E-Lj  ( t )  |  and  e  =  Max  (2.2) 

t  i,j 


Thus  e  is  the  magnitude  of  the  maximum  deviation  expected  in  the  entries  of  A. 
In  the  case  of  'unstructured'  perturbations  only  a  'norm'  of  the  matrix  E,  say 
°rmax^E)  i-s  assumed  to  be  known. 

J _ ^2_uH^s  for  Structured  Perturbation 

For  this  situation,  it  is  shown  in  [13],  [14]  that  the  system  of  (2.1)  is 
stable  if 


e 


1 


< - 

n  amax  tPl 


eP 


(2.3) 
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(2.4) 


where  P  is  the  solution  of  the  Lyapunov  matrix  equation 
T 

AqP  +  PA0  +  2In  =0  (In  is  an  nxn  identity  matrix) 
and  it  is  assumed  that  Ey  t  0  for  all  i,j  =  l,2...n. 

At  this  stage,  it  is  important  to  point  out  that  Patel  and  Toda  [14]  do 
not  formally  exploit  the  distinction  between  the  different  types  of 
perturbations  as  discussed  in  the  Introduction  and  employ  a  single  measure  as 
the  basis  for  an  upper  bound  for  all  these  types  of  perturbations.  In  this 
research,  we  clearly  delineate  this  distinction  and  provide  an  improved  upper 
bound  for  structured  perturbations. 

In  what  follows,  we  present  the  main  mathematical  result  to  get  an  upper 
bound  for  the  perturbation  matrix  E  of  system  (2.1),  assuming  structured 
perturbation. 

Theorem  2.1:  The  system  of  (2.1)  is  stable  if 
1 

£  < - — — — —  -  y  (2.5) 

amax  fpm^n^s 

where  Un  is  an  nxn  matrix  whose  entries  are  unity  i.e.,  Un-j_-j  =  1  for  all  i,j  = 
1 , • . «n  and  P  satisfies  the  Lyapunov  equation  given  by  (2.4)  and  Pm  is  the 
matrix  formed  by  the  modulus  of  the  entries  of  matrix  P  i.e.,  Pmij  =  lpij  I* 

Proof:  Given  in  Appendix  A. 

However,  the  above  expression  for  is  obtained  without  using  the 

'structural'  information  about  the  nominally  stable  matrix  AQ.  Obviously,  one 
can  get  an  even  better  bound  if  the  structural  information  of  the  matrix  AQ  is 
taken  into  consideration.  For  example,  as  indicted  by  equation  (2.2),  if  one 
knows  the  individual  element  bound  e for  all  i,j,  then  this  information 
can  be  incorporated  into  the  calculation  of  the  bound  and  in  such  cases,  the 
bound  ^ y  takes  the  following  form. 
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Theorem  2.2:  The  system  of  (2.1)  is  stable  if 

1 


e  < 


=  u 


o  ( pmu,= )  <, 

max  m  e  s 


Ye 


(2.6a) 


where  Ue  is  an  nxn  matrix  whose  entries  are  such  that 


Jeij 


eij 

e 


(2.6b) 


Thus  Ueij  =  0  if  the  perturbation  in  A0ij  is  known  to  be  zero  (i.e.,  =  0). 

Similarly  Ue^j  =  1  if  one  doesn't  explicitly  know  eij  for  some  i,j  thereby 
accommodating  the  worst  case  situation.  Hence,  it  can  be  seen  that 

0  <  ueij  <  1  (2.7) 

For  example,  if  the  nominally  stable  matrix  is  diagonal  to  start  with  and 
if  it  is  known  that  there  can  be  no  perturbation  possible  in  the  off  diagonal 
elements  then  Ue  is  taken  to  be  In.  Similarly  if  the  nominally  stable  matrix 
AQ  is  such  that 

a1  3 

0  (2.8a) 

a33 

and  it  is  known  that  there  cannot  be  any  perturbations  in  a23,  a3i  and  a32 
then 


a1 1 

a1  2 

a2 1 

a22 

0 

0 

Ue  = 


1  1 
1  0 
0  1 


(2.8b) 


Similarly,  if  one  also  had  additional  information  that  a^  min  <  a^j  <  aj_j  max, 


then  one  can  calculate  |  =  max  (  |  ay  max  -  ay  )|,  |  (a 


I}  min  ”  aij  )  I  )  for 


ail  ifj  and  determine  lAa^j  |max.  Accordingly,  the  matrix  Ue  elements  will  be 
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(2.9) 


Aaij  I  max 

Thus,  the  above  expression  (2.6)  for  Pye  garners  the  structural  information  of 
the  nominally  stable  (and  thereby  the  perturbation)  matrix  in  a  unified  form, 
without  the  need  for  different  expressions  for  different  situations. 

Remark  2.1:  Since  amax  (PmUe)s  <  tfmax  (PmUn)s  it  can  be  seen  that 

PYe  >  yeY  (2.10) 

Remark  2.2  :  The  location  of  the  perturbation  (reflected  by  the  location 
of  nonzero  entries  in  the  Ue  matrix)  clearly  influences  the  perturbation  bound. 
Table  1  illustrates  this  point. 

2.2.2.  Bounds___f o£ JJns  true  tured perturbations 

Now,  we  present  similar  'stability  robustness*  measures  assuming 
unstructured  perturbations,  namely  perturbations  where  only  the  spectral  norm 
bound  for  the  error  matrix  (of  a  given  model  structure)  is  known  [16]. 
Robustness  Measure  of  Patel  and  Toda  [14] 

For  this  case,  the  condition  for  stability  is  given  by 

1 

°"max  (E)  =  ||e||s  < - =  ^PW  (2.11) 

° max ( p ) 

where,  as  before,  P  satisfies  equation  (2.4).  It  is  also  shown  that 

^  PW  ~  | A  (A0 )  | m^n 

if  Aq  is  normal  thereby  avoiding  the  solution  of  Lyapunov  matrix  equation. 
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Table  1 


Elements  of  A  in  which  perturbation  is  assumed, 
o 


Let 


A  = 
o 


-3  -2 

1  0 


li 


ID 


“11 

only 


12 

only 


21 

only 


a. 


22 

only 


all  & 


12 


all  & 


22 


11 


21 


12 

a21 


al2  & 


22 


a21  & 


22 


an' 

ai2' 

*21 


11' 

*12' 


22 


all' 

a2l' 


22 


12' 

*21' 


22 


1  1 
u  u 


i  o 
o  o 


o  1 

o  o 


0  0 
1  0 


0  0 
0  1 


1  1 

0  0 


1  0 

0  1 


0  1 
1  0 


0  1 
0  1 


0  0 
1  1 


1  1 
1  0 


1  1 
0  1 


1  0 
1  1 


0  1 
1  1 


ye 


0.236 


1.657 


1.657 


0.655 


0.396 


1.0 


0.382 


0.40 


0.5 


0.324 


0.3027 


0.397 


0.311 


0.273 


0.256 


Robustness  Measure  ofMLee  [17 ] 

Assuming  that  the  orthogonal  matrix  U  in  the  polar  decomposition  of  A, 

Aq  =  UHR  or  A0  =  HlU 

is  a  stable  matrix,  Lee  gives  a  stability  condition  as  follows: 

^max(S)  <  “CTmin(A0)  cos(0min)  =  ^LW  (2.12) 

where  0min  is  the  smallest  principal  phase  of  A0  measured  counter  clockwise 
from  the  positive  real  axis.  Also  ULW  -  I  X(A0)|min  tf  A0  normal. 

In  the  above  two  measures,  the  measure  Ppw  doesn't  exploit  other  forms  of 
•structural'  information  about  AQ  (except  for  the  case  of  'normal'  matrix) 
whereas  the  measure  ^Lw  imposes  a  restriction  on  A  (which  in  turn  requires 
testing  yet  another  condition).  In  what  follows,  exploiting  other  structural 
information  about  the  nominally  stable  matrix  AQ,  we  present  a  new  stability 
robustness  measure  for  unstructured  perturbations. 

AQ+ Ao 

Theorem  2.3:  Let  the  nominal  matrix  Aq  be  stable  and  let  As  =  -  be 

2 


negative  definite.  Then  the  matrix  Ag  +  Es  is  negative  definite  and  hence 
(AQ+E)  is  stable,  if 

amax^E^  <  °min^As^  =  ^YW  (2.13) 

Proof :  Given  in  Appendix  B. 

Note  that  there  is  no  need  to  solve  the  Lyapunov  equation  to  find  Myw* 

The  following  example  compares  the  bounds  JJpw,  and  yYW  when  the  given 

nominal  matrix  is  stable  and  As  is  negative  definite  (and  Aq  is  not  normal). 

Ex  ample  2.1 


A 


o 


-1  -0.25 

+0.25  -1.2 
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The  eigenvalues  of  AQ  are  -1.1  ±  j  0.229 


MPW  PLW  yYwt=  ^min(As>] 


1.0025  1.0025  1.0 

In  the  above  example  even  though  ^  is  less  than  Upw,  there  is 
considerable  savings  in  the  computation  in  its  determination. 

It  may  thus  be  seen  that,  depending  on  the  structural  information 
available  about  AQ,  one  can  obtain  different  bounds.  For  completeness,  we 
summarize  in  Table  2,  the  available  bounds  for  unstructured  perturbations  and 
the  related  conditions  under  which  they  are  obtained. 

The  main  contribution  of  pertains  to  Case  (3)  in  the  table.  As  noted, 

when  the  nominally  stable  matrix  AQ  is  not  normal  but  As  is  negative  definite, 
one  can  get  ^yw  quite  easily  wheres  ^Pw  and  ^Lw  require  considerable 
computation.  Of  course,  the  bound  ^y^  is  more  conservative  than  because 
CTmin(As)  <  yPW  as  shown  by  Patel,  Toda  [14]  and  Yasuda,  Hirai  [18].  But 
example  2.1  above  clearly  demonstrates  that,  in  some  cases,  ^yw  can  offer 
appreciable  savings  in  the  computation  without  much  sacrifice  in  the  bound. 

It  may  be  helpful  to  mention  here  that  one  necessary  condition  on  a  stable 
matrix  AQ  such  that  As  is  negative  definite  is  that  a^i  <  0  for  all  i  = 

1,2..  .n .  Incidentally,  a  sufficient  condition  is  that  A0  be  normal.  In  fact, 
another  coefficient  condition  can  be  proposed  by  applying  theorems  2.1  and  2.  2 
of  Section  2.2.  Writing 

2As  =  Ao  +  AqT  (2.14) 

and  treating  AT  as  an  'error'  matrix,  we  observe  that  Ag  is  stable  (i.e.,  is 
negative  definite)  if 
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Table  2.  Summary  of  Bounds  for  Unstructured  Perturbations. 


Nominally  stable 


Bounds  on  the  error 


i 


1) 


2) 


Matrix  A 

o 

norm| 

[ E | |  for  stability 

Remarks 

AQsf.A  Symmetric  (normal 
and  negative  definite 

y  - 
PW 

ULW  =  WYW 

Vn(flo> 

Simplest  case  but  restrictive. 
No  need  for  Lyapunov. 

A^  is  negative  definite 

and  A  is  normal 

o 

It  II 

£  S 

ft  X 

p  p 

y  =  1  X  (A)  1 

LW  1  c  1  rran 

a  .  (A  ) 
nun  s 

No  need  for  Lyapunov  in  all 
the  cases. 

A^  is  negative  definite, 

but  A  is  not  normal 

o 

yPW 

1/0  (P) 

max 

Lyapunov  solution  needed. 

PLW  " 

-a  .  (A  )  cos0  . 
min  o  min 

Polar  decomposition  U  to  be 
stable. 

yYW  = 

a  .  (A  ) 
min  s 

Lyapunov  solution  not  needed. 

General  A  (not  normal) 
o 

and  A^  is  not  negative 

definite 

yPW  = 

l/o  (P) 

max 

Lyapunov  solution  needed. 

yLW 

-a  .  ( A)  cos0  . 

min  «  mm 

Polar  decomposition  U  to  be 
stable. 

No  y  possible 
YW 


I  A  T. .  I 

lrto  13  I  max  <• 


amax 


(2.15) 


where  P  satisfies  PAQ  +  AoTP  +  2In  =  0  and  U€ 


I  T 

t^o  ij 


1  ,  2  •  •  .  n  . 


ij  I  a  T .  .  I 

I  ij  I  max 


for  all  i,j  = 


2.2.3  Application  to  an  Ai rcra f t_Con trol  Problem 

We  now  consider  the  same  application  example  as  the  one  considered  by 
Patel,  Toda  &  Sridhar  in  [15].  For  completeness  sake,  we  briefly  reproduce 
here  the  mathematical  model  of  [15] • 

In  [15],  the  system  chosen  is  the  flare  control  of  the  Augmentor  Wing  Jet 
STOL  Research  Aircraft  (AWJSRA).  The  purpose  of  the  flare  control  is  to  make  a 
smooth  transition  from  an  initial  steep  flight  path  angle  of  -7.5°  on  the  glide 
slope  at  an  altitude  of  approximately  65  ft.  to  a  final  smaller  flight  path 
angle  (-1°)  more  appropriate  for  touchdown. 

The  equations  for  the  longitudinal  dynamics  of  the  AWJSRA  at  an  airspeed 
of  110  ft/s  and  flight  path  angle  of  -1°  are  given  by 

x  =  A  x  +  B  u  (2.16) 

where  xT  =  [6  v  6y  69  5q  5h]T 

uT  =  [6e  6n]T 
6 v  =  change  in  airspeed,  ft/s 
6y  =  change  in  flight  path  angle,  deg 
50  =  change  in  pitch  angle,  deg. 

5q  -  change  in  pitch  rate,  deg/s 
6h  =  deviation  from  nominal  altitude,  ft 
6e  =  change  in  elevator  deflection,  deg 
6n  =  change  in  nozzle  angle,  deg 
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-0.0547 

-0.298 

-0.2639 

-0.0031 

0.0^ 

0.16 

-0.4712 

0.4661 

0.0437 

0.0 

A  = 

0.0 

0.0 

0.0 

1  .0 

0.0 

(2.17a) 

0.1751 

0.1236 

-0.1236 

-1  .3 

0.0 

0.0174 

1  .92 

0.0 

0.0 

0.0 

-0.00315 

-0.0943“ 

0.0408 

0.0224 

B  = 

0.0 

0.0 

(2.17b) 

-1 .1200 

-0.08 

_  0.0 

0.0 

The 

open  loop  poles 

of  the  system  are  at 

0.0,  -0.0105 

±  j  0 . 2737 ,  -0.6757 

and  -1.129,  The  control  gain  is  obtained  by  minimizing  [19]  the  performance 
index 


J  ~  J  (xtQx  +  uTRu)dt  (2.18a) 

o 

with  R  =  Diag  [16,  0.5]  and  Q  =  ql5.  (2.18b) 

Applying  the  analysis  of  [15]  and  this  paper,  the  bounds  Uy  (with  Un)  and 
and  their  variation  with  q  are  summarized  in  Table  3  and  Fig.  1. 

Clearly  it  is  seen  that  ^y  is  greater  than  ^p  for  the  values  of  q 
considered  and  the  'optimism'  of  W y  over  increases  as  q  is  increased. 


Table  3.  Variation  of  ^y  and  ^p  with  q 


q 

UP 

yY 

0.1 

0.0055 

0.0061 

0.25 

0.0082 

0.0093 

0.5 

0.0107 

0.0125 

1  .0 

0.01 37 

0.0164 

5 

0.021 3 

0.0272 

10 

0.0240 

0.0322 

50 

0.0305 

0.0420 

102 

0.0323 

0.0451 

104 

0.0364 

0.0530 
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q 


Fig.  1.  Plot  of  |j  vs .  q 
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2 . 3  Reduction  in  Conservatism  by  State  Transformation : 

The  proposed  stability  robustness  measures  presented  in  the  previous 
section  were  basically  derived  using  the  Lyapunov  stability  theorem,  which  is 
known  to  yield  conservative  results.  The  'improvement'  obtained  in  the 
proposed  bounds  is  the  result  of  exploiting  the  'structural'  information  about 
the  perturbation.  Clearly,  another  avenue  available  to  further  reduce  the 
conservatism  is  to  exploit  the  flexibility  available  in  the  construction  of 
Lyapunov  function  used  in  the  analysis.  In  this  section,  a  method  to  further 
reduce  the  conservatism  on  the  element  bounds  (for  structural  perturbation)  is 
proposed  by  using  state  transformation.  This  reduction  in  conservatism  is 
obtained  by  exploiting  the  variance  of  the  'Lyapunov  criterion  conservatism' 
with  respect  to  the  basis  of  the  vector  space  in  which  the  function  is 
constructed.  The  proposed  transformation  technique  seems  to  almost  always 
increase  the  region  of  guaranteed  stability  and  thus  is  found  to  be  useful  in 
many  engineering  applications. 

2.3.1  State  Transformation  and  its  Implications  on  Bounds 

It  may  be  easily  shown  that  the  linear  system  (2.1)  is  stable  (or 
asymptotically  stable)  if  and  only  if  the  system 

A  A  A 

x ( t )  =  A(t)  x(t)  (2.19a) 

where 

x(t)  =  Q-1  x(t),  A(t)  =  Q-1  A(t)Q  (2.19b) 

and  Q  is  a  nonsingular  time  invariant  nxn  matrix,  is  stable  (or  asymptotically 
stable ) . 

Even  though  the  proof  of  this  result  is  quite  straightforward,  it  is  to  be 
emphasized  that  it  is  not  based  on  the  standard  eigenvalue  argument  as  we  are 
dealing  with  time  varying  case. 
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The  implication  of  this  theorem  is,  of  course,  important  in  the  proposed 
analysis  here.  It  means  that,  to  investigate  the  stability  of  a  linear  system 
of  (2.1),  one  can  transform  it,  by  a  linear  map,  to  a  different  coordinate 
frame  and  derive  the  stability  robustness  condition  in  the  new  (transformed) 
coordinates.  Since  the  available  stability  robustness  conditions  are  just 
sufficient  conditions  (in  contrast  to  necessary  and  sufficient  conditions),  it 
is  likely  that  one  can  get  a  less  conservative  bound  in  the  transformed 
coordinates . 

In  sequel,  we  consider  a  diagonal  transformation  matrix  Q.  (The  use  of  a 
general  non-diagonal  transformation  matrix  is  under  study. )  Let 

Q  =  diag [q^ ,q2 , • . .qn]  qi  *  0,  i  =  1,2,...,n  (2.20) 


Then  A(t)  =  Aq  +  E(t)  = 
ai i+ei i 

^1 /q2 ( a21 +e21 (t) ) 

qiAln(an1+<2n1  <t) ) 


q2/qi <a1 2+e1 ) . qn/qi (ain+e1n(t) ) 

a22+e22( t) . qn/q2(a2n+e2n(t)  > 


(2.21  a) 


q2/qn ^ an2+en2 ( ) . (ann+finn(t)  ) 


where  A0  =  Q~^A0Q  and  E(t)  =  Q  ^E(t)Q. 
Correspondingly,  we  get 


eij 


and  e 


Max  * 

.  eij 

h] 


(2.21b) 


(2.22) 


and  as  discussed  before 


U, 


ei] 


-  e  i  j  /  e 


(2.23) 


It  may  be  seen  that  (2.21)  is  similar  in  form  to  the  'weighted  norm' 
matrix  which  has  been  used  successfully  in  frequency  domain  to  reduce  the 
conservatism  of  the  stability  robustness  condition  [20,  21,  22]. 
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In  general ,  stability  robustness  conditions  of  the  type  given  in  (2.6)  may 
be  used  in  two  different  situations;  (a)  given  the  perturbation  element 
bounds,  e-Lj  (of  the  l.h.s.  of  equation  (2.6),  the  condition  (2.6)  can  be  used 
to  check  the  stability  of  the  perturbed  system  (2.1)  as  well  as  check  the 
correponding  conservatism  of  the  condition,  or  (b)  simply  treat  the  condition 
(2.6)  as  specifying  a  bound  on  the  elemental  perturbation  (i.e.,  of 

l.h.s.  of  (2.6)  are  not  explicitly  known  but  Ue  of  r.h.s.  of  (2.6)  is  known? 
recall  that  in  the  absence  of  any  explicit  and  relative  information  on  ,  one 
can  take  Ue-  =  1  thereby  accounting  for  the  worst  case  situation).  This  type 
of  delineation  is  useful  in  arriving  at  a  transformation  appropriate  for  the 
purpose  at  hand.  We  illustrte  these  situations  by  way  of  simple  examples. 

(a)  L.H.S.  is  known,  (checking  stability  and  conservatism)  : 

For  this  case,  in  the  transformed  coordinates,  the  stability  condition 
(2.6)  bcomes 

A  A  A 

£  i j  <  jj  Ue*pj  (2.24a) 

or 


where 


£  <  y 


y  -  l/tfmax  fpm^e^i 


(2.24b) 


(2.25a) 


and  P  satisfies  PA0  +  AqP  +  2In  =  0  (2.25b) 

The  conservatism  of  the  condition  with  respect  to  the  transformation  can 
clearly  be  compared  by  using  the  index,  defined  by 


A  y  -£ 
0  =  — 


*  A  M-e 
B  =  — 


(2.26) 


where  it  can  be  seen  that  3  and  3  >  0  when  the  stability  condition  gets 
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satisfied  and  that 


3  >  3 


(2.27) 


indicates  the  transformation  to  be  effective  in  reducing  the  conservatism  of 
the  condition.  Let  us  illustrate  this  by  an  example. 

Example  2.2:  Let  A0  of  (2.1)  be  given  by 


A 


o  “ 


-3 

1 


-2 

0 


Suppose  that  only  element  an  gets  perturbed  and  that 
e  1 1  =  le 1 1  I  max  =  ^ 


Then  clearly 


e  =  2  and  Ue  = 


1  0  | 
o  oj 

The  r.h.s.  of  equation  (2.6)  gives  y  =  1.657.  Thus 
1 .657  -  2 

8  - - -  -  0.207 

1  .657 


is  negative  indicating  that  the  stability  condition  is  not  satisfied, 
employing  the  transformation  matrix 
Q  =  diag  [1 ,  1000] 

the  transformed  quantities  are 

^  A  A 

£  =  2  and  Ue  -  Ue  and  y  =  3. 

Thus 


-  3-2 

3  = - =  0.333 

3 


However, 


indicating  that  the  stability  condition  is  satisfied  and  hence  the  system  (even 

in  the  original  coordinates)  is  stable.  Thus  the  transformation  reduced  the 
conservatism  of  the  stability  condition. 
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(b)  L .  H .  S  . _ is  not  known _ ( specifying  the  perturbation  bound)  : 

As  before  in  case  (a),  after  the  transformation,  the  stability  condition 
is 


e  <  V  (2.28) 

A  A 

However,  as  a  bound,  y  is  given  on  the  perturbation  e  and  not  on  e  ,  the 
perturbation  in  the  original  coordinates.  Evidently,  in  order  to  examine  the 
usefulness  of  the  transformation  in  getting  an  improved  bound  it  is  necessary 
to  obtain  the  bound  on  e  after  the  transformation.  Let  y*  denote  the  bound  on 
e  after  the  transformation.  Then  we  have  the  following  result: 

Theorem  2.4  :  Given  Q  =  diag  [q-|  ,q2 ,  • . .  , qn]  ,  and  eq.  (2.22)  the  system  of 
(2.1)  is  stable  if 

e  <  P*  (2.29a) 


<5r 


p*  =  p - 


u 


ers 


^s 


where  y  is  given  by  (2.25)  and  r,s  are  such  that  rs 

A  A 

Ue  corresponding  to  Ueij  max  and  Ue  =  Qm-1UeQm. 
Remark  2.3:  Clearly  y*  >  y  indicates  the  reduction 
condition. 

Remark  2_^4:  It  can  also  be  noted  that 
M*  >  V 

if  and  only  if  (proof  is  available) 


(2.29b) 

is  the  specific  entry  in 

in  conservatism  of  the 


(2.30) 


@  >  6 

A 

where  0  and  0  are  given  by  (2.26). 
Example  2.3:  Let  us  again  consider 

-3  -2~ 

1  0 


Ao  = 


(2.31 ) 
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and  let 


U 


e  _ 


1 

1 


1 

0 


Then  we  have  y  =  0.3972. 


By  using  the  transf ormation  Q  =  diag  [1,  1.8]  it  can  be  shown  that 
y*  =  0.473 


indicating  around  20%  improvement  in  the  bound. 

2*3. 2  Application  to  VTOL  Ai rc raf_t  Control 

In  this  section,  we  extend  the  proposed  method  of  reducing  the 
conservatism  by  state  transformation  to  the  VTOL  aircraft  control  example 
discussed  in  [23] .  This  particular  aircraft  control  problem  has  been  used  in 
the  illus trtion  of  several  control  design  methodologies.  For  example  [23] 
suggests  an  adaptive  control  algorithm  for  the  given  range  of  perturbation  while 
[29]  recommends  a  nonlinear  feedback  control  to  achieve  boundedness  which, 
however,  still  cannot  guarantee  robustness  for  the  entire  range  of  perturbation 
considered.  In  the  present  discussion,  we  propose  a  simple  constant  linear 
rate  feedback  to  guarantee  asymptotic  stability  using  the  proposed  perturbation 
bound  analysis  with  state  transformation. 

The  linearized  model  of  the  VTOL  aircraft  in  the  vertical  plane  is 
described  by 

x(t)  =(A0+AA(t))x(t)  +(B+AB(t))u(t)  (2.32) 

The  components  of  the  state  vector  x-*-r4  and  the  control  vector  u*R^  are 
given  by 

X1  +  horizontal  velocity  (knots) 
x2  +  vertical  velocity  (knots) 

X3  pitch  rate  (degrees/sec) 
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x4  "*■  pitch  angle  (degrees) 
ui  'collective '  pitch  control 

u2  ■>  'longitudinal  cyclic'  pitch  control 
In  [23]  ,  it  is  shown  that  significant  changes  take  place  only  in  the 

elements  832,  334  and  b2i  •  The  ranges  of  values  taken  by  these  elements  are 
given  in  [23]  as 


0.0663  <  a32(~0. 3681 )  <  0.5044 

0.122  <  a34 (~ 1 . 422 )  <  2.528  (2.34) 

0.977  <  b21 (~ 3 . 544 )  <  5.1114 

Thus  the  perturbation  bounds  are  asymmetric  with  respect  to  the  nominal. 


value.  In  order  to  take  full  advantage  of  the  'Perturbation  Bound  Analysis' 
presented  in  Section  II,  we  will  'bias'  the  nominal  value  of  a32,  a34  and  b21 
such  that  we  obtain  symmetric  bounds.  Accordingly,  the  nominal  values  of  a32, 
a34  and  b2i  are  a32  =  0.2855,  a34  =  1.3229,  b2i  =  3.04475.  The  full  matrices 
AQ  and  Bq  are  given  by 


-0.0366 

0.0271 

0.0188 

-0.4555 

0.0482 

-1  .01 

0.0024 

-4.0208 

0.1002 

0.2855 

-0.707 

1 .3229 

0 

0 

1 

0 

(2.35a) 


T  _ 


0.4422  3.04475 

0.1761  -7.5922 


-5.52 

4.49 


0 

0 


(2.35b) 


so  that 
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I  2  I  max  -  0.2197 
lAA34|max  ~  1-2031 

I A  Ay  |  =  o  for  all  other  i  and  j  (2.36) 

lAB21 Imax  =  2.06725 

|ABy  |max  =  0  for  all  other  i  and  j 

A  Robust  constant  gain  linear  state  feedback  control  law  is  obtained  as 
follows : 

Step__1j  Assuming  a  quadratic  performance  index 

oo 

J  =  /  (xTx  +  puTu)dt  ,  (2.37) 

o 

where  p  is  the  control  weighting,  the  optimal  feedback  gain  G-j  is  obtained  as 


where 


Gt  =  -  -  bqtk 
p 

K  B0B0TK 

KA0+A0Tk - +  I4  =  0 

p 


(2.38a) 


(2.38b) 


and  the  closed 
A  = 


loop  nominal  matrix 
Aq  +  BqGi 


(2.39) 


is  found  to  be  asymptotically  stable  as  the  conditions  of  complete 
controllability  and  observability  are  satisfied. 

The  perturbed  closed  loop  system  is  then  given  by 

x(t)  =  [  (Aq+BqGi  )  +  (AA+AB  G-|  )  ]  x  (t)  (2.40a) 

=  [  ( A0+B0G-|  )  +  E]  x  (t)  (2.40b) 

where  Em  =  AAm  +  AB^*.  Let  eij  i  ^  . 
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.sjteP__2:  Clearly  in  this  case  the  perturbation  matrix  is  fully  known.  The 

variable  p  is  varied  such  that  the  stability  robustness  index 
A  y  ( G-j  )  -  e  (G-j ) 


is  made  as  large  as  possible  where  jj  is  given  by 

^  “  1/crmax  fpmue^  s  (2.42a) 

and  Pm  satisfies 

PUo+BqGt)  +  (Ao+BoGt)1  P  +2In  =  0  (2.42b) 

(Note  that  both  jj  and  e  are  functions  of  gain  G1  ) . 

The  gain  that  makes  6  of  (2.41)  maximum  for  the  given  perturbations 
(2.36)  is  obtained  as 


-0.467  0.01388  0.539  0.806  *| 

G1  = 

_  0.043  0.3828  -0.1899  -0.5947_J 

and  the  corresponding  value  of  0  is  given  by 

6  =  -  0.21 


(2.43) 


(2.44) 


which  shows  that  the  stability  robustness  condition  is  not  satisfied  with  the 
gain  G1  of  (2.43) . 

Step  3:  Evidently,  one  needs  to  look  for  a  gain  such  that  0  is  positive.  In 
other  words,  one  needs  to  increase  the  bound  p  without  much  increase  in  e  in 
equation  (2.6).  Toward  this  direction,  for  this  particular  example,  we  propose 
to  use  a  gain  G2  such  that 

A 

Aq  +  BqGt  +  BqG2  =  Aq  +  B0G  (G  =  Gi  +  G2)  (2.45) 


is  asymptotically  stable  and 


(2.46) 


=  AAm  +  (i.e.,  GmG2m  =  0) 

One  form  for  G2  that  satisfies  (  )  since  only  B21  is  nonzero)  is 

|  0  0  0  0 

G2  =  | 

1-321  922  923  924 


The  motivation  behind  the  selection  of  a  specific  G2  matrix  as  above  is  to 
make  the  norm  of  the  matrix  Aq+b0G  bigger  with  the  hope  of  decreasing  the  norm 
of  the  Lyapunov  matrix  P  of  (2.42)  which  in  turn  increases  the  bound  y  without 
increasing  e  of  (2.6).  With  this  in  mind,  we  obtain  g2 =  g23  =  0  and  g22  >  0 
and  g24<  0.  For  simplicity,  we  choose  | g 2 2  I  “  | g 24  I • 

To  guarantee  the  asymptotic  stability  of  A0+B0G  of  (2.45)  we  could  think 
of  the  matrix  BQG2  as  a  perturbation  on  the  nominal  matrix  A0+B0G<j  and  apply 
the  perturbation  bound  condition  of  (2.6).  With  this  done,  we  get 


7.6g22  <  1 .0342 

which  makes  g22  =  0.1362  and  g24  =  -0.1362. 

Thus,  we  finally  get 

|  -0.467 
G  ~  G1 +Go  =  I 

[__  0.043 

and 

'“-0.2356 


0.01388  0.539 

0.519  -0.1899 


Ao+BoG  = 


-1 .7021 
2.8732 
0 


0.1 246 
-4.908 
2.539 
0 


0.22377 

3.0859 

-4.5359 

1 


0.806 

-0.731 

-0.2277 

3.98 

-6.408 

0 


(2.47a) 


(2.47b) 


Now  the  computation  of  3  of  (2.41  )  with  this  new  closed  loop  system  gives 

(2.48) 


6new  =  “  0.1722 


-27- 


which  is  still  negative  but  is  an  improvement  over  the  3  obtained  before. 
However,  one  cannot  conclude  the  stability  of  the  new  perturbed  closed  loop 
system  using  this  Bnew. 

?tep_4:  We  now  appiy  State  transformation  for  the  above  nominal  system  using  a 

transformation  matrix 

Q  =  diag  [0.7,  1.6,  1.02,  0.8]  (2.49) 

A 

The  computation  of  3 new  (after  the  transformation)  gives 

A 

3 xi s w  “  0 -029  y  0 

Thus,  the  state  transformation  given  by  the  matrix  Q  of  (2.49)  has  reduced 
the  conservatism  of  the  stability  condition  and  the  gain  matrix  given  by  (2.47) 
guarantees  robust  asymptotic  stability  of  the  system  in  the  entire  range  of 
parameter  perturbations  given  by  (2.34)  (or  (2.36)). 

2 • 4  Reduced  Conservatism  in  Ultimate  Boundedness  Control 

In  this  section,  we  apply  the  concept  of  reducing  conservatism  by  state 
transformation  to  the  problem  of  obtaining  stabilizing  controllers  for 
uncertain  dynamical  systems  using  the  theory  of  ultimate  boundedness  notably 
papers  by  Leitmann,  Barmish  and  their  colleagues  [24,25,26,27,28]. 

The  basic  philosophy  of  ultimate  boundedness  control  is  that  it  guarantees 
the  state  of  the  system  to  enter  and  remain  within  a  prescribed  neighborhood  of 
the  origin  (zero  set  point)  or  even  a  nonzero  set  point.  In  refs.  [24-25] 
controllers  are  proposed  that  guarantee  ultimate  boundedness  for  systems 
satisfying  the  so-called  'matching  conditions'  namely,  certain  preconditions  on 
the  location  and  magnitude  of  the  uncertainty  within  the  system  matrices. 

Under  the  matching  condition  assumption  uncertainties  with  an  arbitrarily  large 
prescribed  bounds  can  be  tolerated.  Barmish  [26]  extended  this  theory  to  more 
general  situations  by  introducing  the  concept  of  'Measure  of  Mismatch’  (and 
'Mismatch  Threshold'). 
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In  this  framework ,  a  certain  'decomposition*  is  performed  on  the  system 
dynamics  which  yields  a  'matched  portion' .  It  is  shown  that  effective  control 
is  possible  as  long  as  the  'measure  of  mismatch*  does  not  exceed  some  critical 
'mismatch  threshold'.  Even  though  this  'mismatched  uncertainty'  increases  the 
applicability  of  the  theory  of  ultimate  boundedness  to  more  general  forms  of 
uncertainty,  it  turns  out  that  the  resulting  bounds  are  very  small  and  thus 
conservative.  The  reason  for  this  conservatism  is  twofold.  Firstly,  the 
decomposition  of  the  system  into  matched  and  mismatched  parts  is  nonunique;  the 
'better'  the  decomposition,  'larger'  the  size  of  the  uncertainty  that  can  be 
tolerated.  Secondly,  the  'mismatch  threshold'  is  calculated  based  on  'Lyapunov 
stability  theory'  which  is  always  known  to  be  conservative  in  predicting 
stability,  as  discussed  in  the  previous  section. 

The  main  focus  of  this  section  is  to  address  these  aspects  of  i) 
consevatism  in  the  bounds  obtained  in  time  domain  using  Lyapunov  theory  and  ii) 
the  nonuniqueness  of  the  decomposition.  Specifically,  a  diagonal  similarity 
transformation  on  the  system  is  proposed  which  serves  to  achieve  either  a  less 
conservative  bound  on  the  mismatched  perturbation  or  a  'better'  decomposition 
of  the  matched  and  mismatched  parts  or  both.  The  proposed  extension  can  be 
applied  to  both  zero  set  point  as  well  as  nonzero  set  point  control  [29] . 

2.4.1  System  D e sc ription  and  Tran sformation  on__Bou ndedness  Set 

As  in  refs.  [26,29],  consider  the  uncertain  dynamic  system  described  by 
the  state  equation 

x ( t )  =  [ A+A  A(r(t) )]x(t)  +  [ B+A  B ( r ( t ) ) ]u(t) +Cv ( t )  (2.50) 

where  x(t)-^Rn  is  the  state,  u(t)+Rm  is  the  control  and  v(t)eRs  is  the 
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disturbance  input.  A,B,C  are  nominal  system  matrices  of  appropriate  dimensions 
and  AA  and  AB  are  the  uncertainty  matrices  depending  continuously  on  the 
uncertain  parameter  vector  r(t)*Rp.  Also 
r  (  •  )  >  fi  rP 
v  ( • )  V  +  RS 

are  Lebesque  measurable,  where  ft  and  V  are  prescribed  compact  subsets  of 
appropriate  spaces.  It  is  also  assumed  that  the  pair  (A,B)  is  completely 
controllable.  Let  us  also  define 

z(t)  =  x(t)  -X*,  x(tQ)  =  xQ 

A  (2.51) 

Z<V  ~  2 o  =  xo  ”  x* 

where  x*  is  the  set  point  of  the  state.  We  denote  the  solution  of  (2.50)  by 
x(t,xQ,t0).  Let  ||(*)|  |  denote  the  Euclidean  norm  for  a  vector  and  spectral 
norm  [=Amax ( • )T( • ) ]  for  a  matrix  (•). 

We  now  restate  the  definitions  used  in  [26]  for  continuity  and  better 
exposition. 

Definition  1 :  Given  a  solution  x(t,x0,t0),  [t0, t^ ]  of  (2.50)  and  a  set  point 

x*+Rn,  we  say  that  the  solution  is  uniformly  bounded  if  there  is  a  positive 
constant  d(zQ)  <  « r  possibly  dependent  on  z0,  but  not  on  tQ,  such  that 
|  |x( t,x0,t0)-x* I  I  <  d(zQ)  for  all  t+[t0,til. 

Definition  2:  Given  a  solution  x(t,x0,t0),  t>(t,«>)  of  (2.50)  and  a  set  point 
x**Rn,  we  say  that  the  solution  is  uniformly  ultimately  bounded  with  respect 

A  A 

to  a  set  S  if  there  is  a  nonnegative  constant  T(z0fS)  <  00 ,  possibly  dependent 

A  A  A 

on  zQ  and  S  but  not  on  tQ/  such  that  x(t,x0,t0)  +S  for  all  t  >  t0+T  (zQ,S)  and 

A 

that  x**S. 
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We  now  present  the  following  theorems  which  are  crucial  in  -justifying  the 
use  of  a  transformation  for  later  analysis. 

^®2.r®5L  “  fx-j  ( t)  ,X2( t) ,  .  .  .  ,xn(t)  ]T  is  uniformly  bounded  to  the  set 

point  x*  =  [x*t ,x*2, . . • ,x*n]T  if  and  only  if 
x(t)  =  [x1 (t) ,X2(t) , . . . ,xn(t) ]T 

=  [k-|Xi  (t)  ,k2x2(t) , .  . .  ,knxn(t)  ]T  (2.52) 

kj  =  constant  *  0  V  i=l,...n. 
is  uniformly  bounded  to  the  set  point 

x*  =  [k^ x* i ,k2X*2 » • • • ,knx*n ] T 

Theorem  2.6:  x(t)  is  uniformly  ultimately  bounded  with  respect  to  a  set 
~  A 

S  =  (x-»-Rn:  (x(t)-x*)T  P  (x  ( t)  -x*  )  <  K)  (2.53) 

where  P  is  a  symmetric  positive  definite  matrix,  if  and  only  if  x(t)  of  (2.52) 

is  uniformly  ultimately  bounded  with  respect  to  a  set 

~  A  _  ... 

S  =  ( x->-  Rn :  (x(t)-x*)T  P  (x  ( t )  —x*  )  <  K)  (2.54) 

These  two  theorems  aid  us  now  to  establish  the  transformation 

x  =  Qx  (2.55) 

where  Q  =  diag  [...qi...],  qi  =  1/ki  V  i  =  1 . . .n  to  be  a  similarity 
transformation.  With  this  transformation  the  system  of  (2.50)  is  transformed 
to 


x(t)  =  (A+AA(r(t) )x(t)  +  ( B+A  B ( r ( t ) ) u ( t ) +Cv ( t ) 

where 

A  =  Q-1AQ,  A A  =  Q-1  AAQ 

B  =  Q-1B,  AB  =  Q-1  AB  and  C  =  Q_1C 


(2.56) 


(2.57) 
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We  now  recall  from  theorems  2.5  and  2.6  that  if  the  solution  of  (2.50)  is 
uniformly  bounded  and  uniformly  ultimately  bounded  with  respect  to  a  specific 
set,  then  so  is  the  solution  of  (2.56).  We  also  note  from  [26]  that  the  main 
condition  for  a  solution  of  (2.50)  to  be  bounded  is  the  so  called  mismatch 
constraint  condition 


M  <  M*  (2.58) 

where  M  is  the  'measure  of  mismatch1  and  M*  is  the  ’mismatch  threshold'.  The 
motivation  behind  employing  the  transformation  (2.55)  is  to  basically  transform 
the  above  mismatch  condition  with  the  hope  that  it  would  yield  an  improvement 
over  the  original  condition  (2.58)  in  the  following  sense: 

Define  g,  the  'Mismatch  Conservatism  Index'  (MCI)  in  the  original  system 


as 


A 

6  «  (M*  -  M)/M* 


and  (3,  the  MCI  for  the  transformed  system 


(2.59) 


_  a  _  r  __ 

P  -  (M*  -  M)/M*  (2.60) 

Then  if  8  >  3  ,  we  have  a  less  conservative  mismatch  constraint  condition. 


2.4.2  Control  Strategy  and  Mismatch  Condition  Under  Transformation 

The  nonlinear  control  law  obtained  for  the  original  system  [26]  (for 

simplicity,  taking  the  set  point  x*  to  be  zero)  is  given  by 

u(t)  =  Gx  +  p ( x )  (2.61  ) 

where  the  constant  gain  matrix  G  is  the  linear  feedback  gain  that  makes 
A 

Ac  =  A  +  BG  (2.62) 
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stable.  Following  the  lines  of  [26] ,  the  decomposition  of  the  uncertainty 


matrices  into  matched  and  mismatched  portions  is  given  by 

AA(r(t) )  =  AAm(r(t) )  +  AA(r(t))  =  BD(r(t))  +  AA(r(t)) 

AB ( r ( t ) )  =  ABm(r(t))  +  AB(r(t))  =  BE(r(t))  +  AB(r(t))  (2,63) 

C  =  Cm  +  C  =  BF  +  C 
Provided  that 

E  =  [1  -  Max  | |E(r(t) ) | | ]  >0 
r+R 

we  define 

A  2  =  E"1  [Max|  |D(r(t))|  |  +  Max|  | E (r ( t ) )G |  |] 
r-*R  r*R 


and  furthermore 

M  =  Max | |AA(r(t)| |  +  Max| | AB(r ( t) )G | | 
r+R  r*R 

+  A  2  MaX I  | AB(r ( t ) ) |  | 
r+R 

Then,  the  mismatch  constraint  condition  for  the  original  system  (2,50)  is 

'V/ 

M  <  M* 

where  M*  =  1/(2Amax(p^  (2.64a) 

and  P  is  the  symmetric  positive  definite  matrix  obtained  as  the  solution  of  the 
Lyapunov  equation 

PAC  +  AcTp  +  In  =  0  (2.64b) 

It  may  be  noted  that  M*  as  given  by  eg.  (2.64)  is  exactly  the  same 
quantity  obtained  by  Patel,  Toda  in  [14]  as  the  bound  for  the  nonlinear  (as 
well  as  linear  unstructured)  perturbation  (on  an  asymptotically  stable  linear 
system)  that  guarantees  the  stability  of  the  perturbed  system. 

Now  applying  the  suggested  transformation  (2.55)  on  (2.61)  and  it  gives 

u  =  G  x(t)  +  p(x)  (2.65) 
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and 


x(t)  =  (A+AA(r(t)  )x(t)  +  { B+A B(r(t))(G  x(t)+p(x))  +  Cv(t)  (2.66a) 

where  G~  =  GQ  (2.66b) 

and  p(x)  =  p(Q  x)  (2.66c) 

Let  the  decomposition  be  given  by 

AA(r(t) )  =  AA^rCt))  +  AA(r(t))  =  B  D  (r(t))  +  AA(t(t)) 

AB(r  ( t) )  =  A  Bjh  ( r  ( t ) )  +  AB(r(t))  =  B  E(r(t))  +  AB(r(t))  (2.67) 

C  =  Cm  +  C  =  B  F  +  C 
Provided  that 


E  =  (1  _  max  | |E(r(t) ) | | )  >  0 
r+R 


we  define 


A2  =  E"1  { Max | | D ( r ( t ) | |  +  Max| |E(r(t) )g| | ) 
r+R  r+R 


and  furthermore 


M  =  Max | | AA ( r  ( t ) | |  +  Max | | AB(r ( t) )g| |  +  X2  Max  ||AB(r(t))|| 
r*R  r+R  r^R 

and  the  mismatch  constraint  condition  for  this  transformed  system  is 


M  <  M*  (2.68) 

where 

M*  -  1/ (2Amax(P) )  (2.69a) 

where  P  satisfies 

-  --T  _ 

P  ^  +  AcP  +  In  =  0  (2.69b) 

and 


-34- 


Ac  =  A  +  B  G 


(2.69c) 


Note  that  Ac  =  Q^AqQ  is  stable  whenever  Ac  is  stable. 

Now ,  it  may  be  observed  from  the  above  discussion  that  there  are  two 
avenues  available  for  reducing  the  conservatism  of  the  mismatch  condition  by 
transformation.  One  is  to  get  a  ’better1  decomposition  through  eq.  (2.67)  than 
the  original  decomposition  of  (2.63)  in  the  sense  that  bigger  portions  of  the 
uncertainties  may  be  lumped  into  the  matched  part  in  the  transformed  equations. 
The  other  avenue  available  is  to  get  an  improved  bound  M*  (in  the  transformed 
system)  since  conservatism  of  the  Lyapunov  condition  is  not  invariant  under 
transformation.  However,  it  is  interesting  to  note  that  by  the  specific 
transformation  suggested  in  (2.55),  we  may  be  able  to  achieve  both  (better 
decomposition  as  well  as  better  bound)  whereas  in  some  situations  we  may  achieve 
either  one  of  them.  In  some  other  situations  whatever  gained  in  one  avenue  may 
be  lost  in  the  other.  Thus,  in  order  to  assess  the  final  effectiveness  of  any 
proposed  transformation,  one  has  to  compare  3  and  3  (the  Mismatch  Conservatism 
Index)  defined  earlier. 

Once  the  mismatch  constraint  condition  is  analyzed  in  the  transformed 
system  for  possible  reduction  in  conservatism,  it  is  desirable  to  express  the 
control  law  and  the  ultimate  boundedness  measures  like  the  radius  of  the 
largest  closed  ball  wholly  contained  in  the  set  (denoted  by  n  in  [26]  and  nmax 
in  [29] )  in  the  original  coordinates.  This  aspect  is  discussed  in  the  next 
section. 

2.4.3  Control  Law  and  Ultimate  Boundedness  Measures  w.r.t.  The  Original  System 

It  may  be  observed  that  the  mismatch  constraint  condition  is  related  to  the 
decomposition  of  the  uncertainties  and  thus  is  related  to  the  nonlinear  control 
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law  as  well.  Since  the  decomposition  employed  in  the  transformed  system  to  get 


a  ’better’  mismatch  condition  is  generally  different  from  the  one  obtained  by 
transformation  on  the  original  decomposition,  i.e.,  D,  E  and  F  of  expression 
(2.67)  are  not  necessrily  similar  to  D,  E  and  F  of  expression  (2.63)  before  the 
transformation  Q  (or  Q“1 ) ,  the  parameters  of  the  control  law  in  the  original 
coordinates  obtained  by  inverse  transformation  could  be  different  from  those  of 
the  control  law  in  the  original  system  before  transformation. 

Control  Law 

Suppose  the  control  law  in  the  transformed  system  is  found  as 


(2.70) 


u  =  G*  x  +  p(x) 


then,  by  inverse  transformation,  the  control  law  in  the  original  coordinates, 
is 


U  =  Gnx  +  Pn(x) 


(2.71a) 


where 


Gn  =  G*  Q"1  ,  Pn(x>  =  P(Q-1  x)  (2 

and  the  matrix  A+BG  =  Q( A  +  B  G*)Q-^  is  stable  because  A  +BG*  is  stable  by 


(2.71b) 


design.  The  nonlinear  function  pn(x)  can  be  written  as  [26] 


I |b  pq-1x| I 


T 


p(Q-1x),  if  | |b1PQ-1x| |  >  e 


(2.72a) 


<  £ 


(2.72b) 


where  p(Q_1  x)  =  A-j  +  A2M2  1xll 


and 
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*1  =  [1  -  max  |  |E(r(t)  |  |  l*1  ">ax  |  |Fv|  | 
r+R  v-*V 

A2  =  [1  -  max  | |E(r(t) | | l"1  (max  | | D ( r ( t ) | |  +  max  | |E(r (t) )g| | ) 
r+R  r>R  r*R 

U_1  timate  Boundedne ss  Measures  (U.B,M.  ) 

One  rough  measure  of  the  ultimate  boundedness  set,  used  for  comparisons, 
is  the  radius  of  the  largest  closed  ball  wholly  contained  within  the  ultimate 
boundedness  set  S.  This  is  denoted  by  n  in  [29]  and  nmax  in  Barmish  [26].  Let 
us  denote  it  by  n •  Recall  that  if  the  solution  of  (2.50)  is  uniformly 
ultimately  bounded  with  respect  to  a  set  S(K) 

A 

S ( K )  =  { Rn :  (x(t)-x*)T  P ( x ( t ) — x * )  <  K}  (2.73) 

then  n  and  K  are  related  by  [26] 

K  “  Amax(P)R^  (2.74) 

From  the  theorems  2.5  and  2.6  we  can  observe  that  if  the  solution  of  the 
transformed  system  (2.56),  x(t),  is  uniformly  ultimately  bounded  with  respect 
to  a  set 

S(K)  =  {x+Rn:  (x(t)-x*)TP(x(t)-x*)  <  K}  (2.75) 

where  K  =  Amax  (P)n^  (and  ri  is  the  U.B.M.  in  the  transformed  system)  then  the 
solution  of  the  original  system  x(t)  =  Qx(t)  is  uniformly  ultimately  bounded 
with  respect  to  a  set 

At 

S(Kn)  =  { x+  Rn :  (x(t)-x*)TP(x(t)-x*)  <  Kn  (2.76a) 

where 

Kn  "  Amax )  Amin )  (max  |  qi  |  ) ^  K  (2.76b) 

i 

Note  that  in  (2.76)  the  weighting  matrix  in  the  boundedness  set  with 
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original  coordinates  is  P,  instead  of  P,  the  weighting  matrix  in  the  set  of  the 
original  system  before  transformation.  In  fact,  it  is  one  viewpoint  that 
accounts  for  why  we  can  expect  to  get  a  better  mismatch  constraint  condition. 

From  eqs.  (2.74),  (2.76),  we  observe  that  the  radius  of  the  largest  closed 

ball  wholly  contained  within  the  ultimate  boundedness  set  S(Kn),  denoted  by  nn, 
can  now  be  related  to  Kn  as 

Kn  “  (2.77) 

Thus,  from  (2.75),  (2.76b),  we  get 

nn  =  Umax(p)Amin(p)  ] 1//2  (maxlqj  )n  (2.78) 

i 

Hence,  once  rj  (in  the  transformed  system)  is  known,  one  can  get  r\n*  Of 
course  from  stabilization  and  regulation  (performance)  point  of  view,  it  is 
desirable  to  have  rj  (and  nn)  as  small  as  possible. 

Since  the  configuration  of  the  boundedness  ellipsoid,  characterized  by  P 
of  the  boundedness  set  S  of  (2.75)  and  (2.76),  changes  from  one  coordinate 
frame  to  another,  besides  r\ ,  it  is  appropriate  to  consider  another  parameter 
(say,  v)  as  one  more  rough  measure  on  the  size  of  S.  Clearly,  the  radius  of 
the  smallest  closed  ball  wholly  enclosing  the  ultimate  boundedness  set  would 


serve  the  purpose  (complementing  r\ )  .  Accordingly,  we  define  v 

V  =  t^max ) /^min ^ ^  ^  H  (2.79) 

as  the  radius  of  the  smallest  closed  ball  wholly  enclosing  the  ultimate 
boundedness  set.  Based  on  this  definition  and  (2.75),  (2.76)  we  also  have 


V  »  Umax(p)/Amin(p)]1/2  n  (2.80) 

and 

vn  =  t^max  5Amin^^  hn  (2.81) 
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Roughly  speaking,  the  parameters  n  and  v  serve  as  the  measures  of  minor 
axis  and  major  axis,  respectfully,  in  a  two-dimensional  (ellipse)  setting. 


Origina lCoordinates 
~xft> 

Control  Law 

uo  =  Gx+ p(x) 

G  maximizes  M* 
M.C.I.:  g 
U.  B  .M.  :  r]  ,v 


Transformed  Coordinates  x(t) 

x  ( t )  =  Q_1  x  ( t ) 

Control  Law 

u  =  G*x  +  p(x) 

G*  maximizes  M*(G=GQ*G*) 

M.C.I .  :  {f  >  8 

U.  B .  M.  :  r]  ,  v 


Original  Coordinates 
x(t) 

Control^  Law 
u  =  Gnx  +  Pn(x) 

Gn  =  G*  Q-1 
U.B.M.:  r)n,  vn 


Summary  of  the  Algo r i thm 


1 1 lu s t r a t i ve  Example 


Example  2.4:  Consider 


x-,(t)  =  r-|  ( t)  x-|(t)  +  [2+r2(t) ]x2(t)  +  2v(t) 

x2(t)  =  -  1/2  x*j  ( t )  +  u ( t )  +  v ( t ) 


(2.82) 
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having  uncertainties  bounded  as  follows: 

r(t)  =  {  [ r -|  r 2 1  r  t^+T2^  <  r^  =  constant  >  0 

v(t)  +V  =  {v:  | v |  <  pr  =  constant  >  0 

We  define  the  nominal  pair 


"  0 

2 

"  0~’ 

(A,B) 

-1/2 

0_ 

r 

1 

and  decompose  the  uncertain  portion  as  follows: 


r1 

r2 

"  0 

2 

AA(r )  =  AA(r )  = 

0 

0 

;  C  =  Cm  +  C  = 

1 

+ 

_0 

Let  G  =  [0  g2]  with  g2  <  0. 

* 

It  turns  out  that  max  =  0.22,  with  g2  -  -2.25. 


Now  M  =  max  |  |  A  A  ( r )  |  |  =  max  /(r^  +  r2^ )  =  pr 

^l,r2  r1#r2 


The  mismatch  constraint  condition  is  therefore 

pr  <  M*(g2)r  i.e.,  pr  <  0.22  and  g2  =  -2.25 
We  now  calculate  r]  and  v  (pv=e=1  .0)  and  get 


Pr 

0 

0.05 

0.1 

0.15 

0.2 

T1 

9.14 

11.82 

16.72 

28.6 

99.9 

V 

27.72 

35.86 

50.73 

86.77 

303.1 

Now  let  us  employ  the  transformation 


x  =  Qx  = 


2  0 
0  1 


x 


The  nominal  pair  will  then  be  (A,B)  = 
is 


*  0 

1  ’ 

o' 

_-1 

1_ 

r 

_  1 

and  the  decomposition 
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1 


r1 

-  r2 

- 

0 

1 

A A( r )  = 

2 

=  AA(r ) ;  C  =  Cm  +  c  = 

+ 

.  Now  G  =  [0 

0 

0  „ 

1 

JO 

and 

Max  _  __ 

_  M*  (92)  =  0.3  with  <?2  =  -  1-57 


_  Max  _ _ _  _ 

The  measure  of  mismatch  is  M  =  /  r-^z  +  ( 1/4 ) and  thus 

r1  ,  r2 

M  <  M  and  M*  <  M*  (f  >  g 

Hence,  we  can  conclude  that  the  ultimate  boundedness  of  (2.82  )  is  assured 
if  pr  <  0.3  instead  of  pr  <  0.22  as  given  by  the  analysis  before 
transformation. 

Remark  2.5:  It  is  interesting  to  note  that  the  magnitude  of  the  linear  control 
gain  of  the  final  control  law  (=1.57)  is  not  only  less  than  the  control  gain 
originally  designed  (=2.25),  but  also  tolerates  larger  uncerainties  as 
evidenced  by  bigger  M*.  Thus  it  brings  out  the  point  that  higher  control 
effort  doesn't  necessarily  mean  better  robustness  (where  tolerable  uncertainty 
is  the  measure  of  robustness).  The  calculation  of  rin  and  vn  gives  (py=e"=1  .0) 


Pr 

0 

0.05 

0.1 

0.15 

0.2 

0.25 

0.29 

^n 

14.2 

17 

21 

28 

41 

82 

399 

vn 

29 

35 

43 

57 

85 

169 

820 

Remark  2.6:  Comparing  the  above  tables  and  Figs.  2  &  3,  it  may  be  noted  gn  is 
higher  than  n  (recall  that  both  gn  and  n  are  U.B.M's  in  the  original 
coordinates,  gn  being  the  U.B.M.  after  the  transformation  and  g ,  before  the 
transformation)  in  the  range  0  <  pr  <  0.15,  thereby  indicating  deterioration  in 
the  performance.  But  it  is  to  be  kept  in  mind  that  we  in  turn  gained  in 
stability  in  the  sense  that  the  transformation  allowed  us  to  get  a  'better' 
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mismatch  constraint  condition  thereby  assuring  stability  and  ultimate 


boundedness  for  a  bigger  range  of  uncertainty  (0.2  <  pr  <  0.3).  That  is  n 
doesn't  exist  for  this  range  while  gn  exists.  Thus  the  typical  trade  off 
between  stability  and  performance  is  brought  out  by  this  example.  In  fact,  in 
this  example  there  is  improvement  in  performance  even  for  the  range 
0.15  <  pr  <  0.3. 

Remark  2.7:  Since  in  this  example  M  <  M  and  M*  <  M* ,  we  achieved  reduction  in 
conservatism  from  both  better  decomposition  as  well  as  better  bound  points  of 
view.  Also  note  that  the  particular  transformation  here  may  not  be  the  best  we 
can  use. 

Remark  2.8:  It  may  be  noted  from  this  example  (and  also  from  [26] )  that  the 
gain  which  maximizes  the  mismatch  threshold  M*  also  gives  the  smallest  n  for 
the  same  M  (i.e.,  pr) . 


FIG  2 


nr  ^ 


ROLK  R) 


Figures  2  and  3:  U.B.M.  variation  with  p  . 
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2.5  Application  of  Structured  Bound  to  Stability  of_  Interval  Matrices 

Another  important  contribution  of  the  bound  developed  for  structured 
uncertainty  in  this  research  is  its  applicability  to  the  analysis  of  ’Interval 
Matrices’,  on  which  there  is  interesting  literature  available  [30]-[37]. 

An  interval  matrix  Aj  is  a  real  matrix  in  which  all  elements  are  known 
only  to  belong  to  a  specific  closed  interval  [30] .  An  nxn  interval  matrix  Aj 
is  a  set  of  real  matrices  defined  as  follows: 

Aj  =  [B,C]  =  {A  =  [ay  ]  :  ay+[by,  cy  ]  i,  j  =1  , 2 ,  •  .  .n} 

where 

B  =  [by  ]  ,  C  =  [cy  ]  and  by  <  ay  <  cy  (2.83b) 

It  is  important  to  note  that  there  is  considerable  amount  of  literature 
available  on  the  stability  of  interval  polynomials  starting  from  Kharitonov 
[33],  Bialas  [31],  Bose  [32],  Daoyi  [34].  Bialas  [31]  extends  the  results  of 
the  stability  of  interval  polynomials  to  the  case  of  interval  matrices  and 
presents  necessary  and  sufficient  conditions  for  the  stability  of  interval 
matrices.  But  as  the  counterexamples  given  by  Karl,  Greschak  and  Verghese  [36] 
and  Barmish  and  Hollot  [37]  indicate,  considerable  care  has  to  be  exercised  in 
extending  the  results  of  interval  polynomials  to  the  case  of  interval  matrices. 
Heinen  [35]  and  Daoyi  [34]  propose  simple  sufficient  conditions  for  testing  the 
stability  of  interval  matrices.  While  these  methods  are  indeed  simple  to  use, 
their  application  is  limited  to  the  case  where  ’end  points’  of  C,  namely  of 

equation  (2.83)  are  negative.  In  this  section  we  present  a  different 
sufficient  condition  based  on  the  bound  P Y  of  (2.6)  which  does  not  impose  the 
above  requirements. 
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From  the  results  of  Section  2.2,  it  can  be  easily  concluded  that  given  an 


asymptotically  stable  matrix  F  the  interval  matrix 


[F  -  eUe,  f  +  eUe]  (2.84) 

is  stable  if  e  satisfies  the  condition 

1 

e  < -  (2.85a) 

amax (pmue ) s 

A 

where  Ue  =  [u  =  eij /el  and  P  satisfies  the  Lyapunov  equation  (2.85b) 

eij 

PTP  +  pF  +  2In  =  0  (2.85c) 


In  the  following  section  we  extend  the  above  result  to  establish  the 
stability  of  the  interval  matrix  Aj  of  (2.83). 

2.5.1  Extension  to  Stability  of  Interval  Matrices 

Let  us,  as  before,  denote  the  matrices  A,  B,  C,  etc.  as  follows 

&  ~  ]  ,  B  =  [bj_j  ]  ,  C  —  [cj^j]  ...  i ,  j  =  1 ,2, . .  .n  (2.86) 

It  has  been  shown  by  Bialas  [31]  that  the  interval  matrix  Aj  =  [B,C]  is 
stable  only  if  (necessary  condition)  the  matrices  B  and  C  are  stable.  Thus  we 
start  by  assuming  the  matrices  B  and  C  to  be  stable.  Then  compute  the 
deviations 

A 

dij  =  cij  “  bij  for  a11  if  j  =  1  , 2 ,  .  .  .n  (2.87) 

where  dy  is  seen  to  be  always  nonnegative. 

Denote  D  =  [dy  ]  (2.88) 

A  and  d  =  dij  max  (2.89) 

Form  u  =  dij /d  where  Ue  =  [u  ] .  (2.90) 

eij  eij 

We  then  have  the  following  theorems. 
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Theorem  2.7:  The  interval  matrix  Aj  is  stable  if  the  matrix  B  is  stable  and 
if 


d  < 


1 


amax (p1 mue 


YB 


(2.91 ) 


where 


BtP1  +  +  2In  =  0  and  P 

and  Ue  is  given  by  (2.90). 

Theorem  2.8:  The  interval  matrix  Aj 
if 


lm  is  the  matrix  such  that  P  =  |p 

Imij  1  i  j 

is  stable  if  the  matric  C  is  stable  and 


1 

d  < - 

amax^p2mUe^s 


YC 


(2.92) 


where 

Ctp2  +  P 2 c  +  2In  ”  0  and  p2m  ds  t^ie  matri*  such  that  P  =|p2iii  and 

2mij  *  ' 

U0  is  given  by  (2.90). 

The  above  theorems  clearly  do  not  require  the  diagonal  elements  cii  to  be 
negative  in  contrast  to  the  case  of  Heinen  [35]  and  Daoyi  [34] . 

In  fact,  using  the  above  sufficient  condition  the  interval  over  which  the 
matrix  Ax  is  stable,  can  be  extended  as  illustrated  by  examples  that  follow. 

Finally,  since  both  B  and  C  are  stable  one  can  compute  both  UyB  and  PYC 
and  select  the  one  that  satisfies  either  of  the  conditions  (2.91)  or  (2.92).  If 
both  satisfy  the  condition  ((2.91)  or  (2.92)),  then  one  can  select  the  maximum 
of  the  two  so  that  the  interval  range  can  further  be  expanded, 
matrix  Aj  is  stable,  can  be  extended  as  illustrated  by  examples  that  follow. 
Example  2.5:  Let  us  consider  the  same  example  worked  out  by  Bialas  [31]  and 


-5  1 

-3  2 

Aj  =  [B,C]  where  B  = 

;  C  = 

VO1 

1 

1 

-  5  “4- 
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For  this  problem,  the  condition  of  Bialas  required  extensive  computation  while 
that  of  Heinen  was  not  satisfied.  Let  us  now  apply  the  proposed  technique  to 
this  example. 

Clearly ,  the  matrix  B  is  stable.  Following  equations  (2.88),  (2.89)  and 
(2.90) ,  the  matrices  D  and  are  given  by 


1  2 

i 

1 

i 

i 

in 

. 

o 

D  =  | 

and  thus  d  =  2  and  U  = 

p 

Li 

2__ 

0.5 

1 

Solving  the  Lyapunov  equation 
“0.27273 

P1  = 

0.0909 


(2.91 ) ,  we  get 

0.0909  “| 
0.18182 


and 

UyB  =  2.0736 

Since, 

d  =  2  <  ^ yg  —  2.0736 

we  conclude  that  the  above  interval  matrix  Aj  =  [B,C]  is  stable. 

In  fact,  a  'larger1  interval  matrix  that  is  stable  can  be  obtained  for 
this  example,  as 


A'  =  [B' ,C]  where  B'  = 

l 

0 

;  C  = 

-3 

2 

I 

L  3 

-8__ 

_  5 

-4 

where  Bf  =  B-dUe  and  C  =  B+dUe. 
Example  2.6:  Consider 


-3.2 

-2.2 

|  -3  -2 

Aj  =  [B,C]  where  B  = 

__  0.8 

i 

o 

i 

,  C  =  | 

L 1  0 
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Note  that  c22  =  0  and  thus  none  of  the  above  methods  could  be  applied  to 


this  example.  Applying  the  proposed  technique,  we  have 


D 


0.2 


0.2 


0.2 


0.2 


d  =  0.2  and  thus  Ue 


1 

1 


1 


Noting  that  C  is  stable  and  solving  the  Lyapunov  equation  (2.92),  we  get 


P 


2 


0.5  0.5 

0.5  2.5 


and 


^  Yq  =  0.236 

Since 

d  =  0.2  <  ^YC  =  0.236 


the  given  interval  matrix  is  stable.  A  stable  interval  matrix  with  expanded 
range  can  be  formed  as 


-3.2  -2.2' 

'co 

• 

T 

co 

• 

CsJ 

1 ' 

[B,C ' ]  where  B  = 

;  p*  = 

cJ 

• 

o 

1 

CO 

. 

o 

1 

_  1 .2  0.2 

where 


B  =  C  -  dUe  and  C*  =  C  +  dUe 
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III.  SYNTHESIS  OF  CONTROLLERS  FOR  ROBUST  STABILITY 


The  foregoing  discussion  in  Section  II  is  basically  concerned  with  the 
analysis  of  stability  robustness  for  linear  systems.  No  effort  was  made  to 
synthesize  a  controller  to  achieve  stability  robustness  except  in  2.3.2 
where  the  gain  was  determined  in  an  ad  hoc  way  for  a  specific  problem.  In 
this  section,  we  address  this  design  aspect  from  a  systematic  algorithm  point 
of  view.  The  philosophy  behind  the  proposed  procedure  is  to  make  use  of 
the  perturbation  bounds  developed  in  the  previous  section  in  a  design  formu¬ 
lation  and  give  an  algorithm  to  synthesize  controllers  for  robust  stability. 
Towards  this  direction,  a  quantitative  measure  called  'stability  robustness 
index  is  introduced  and  based  on  this  index  a  design  algorithm  is  presented 
by  which  one  can  pick  a  controller  that  possesses  good  stability 
robustness  property.  The  algorithm,  for  given  size  of  perturbation  can  be 
used  to  select  the  range  of  control  gain  for  which  the  system  is  stability 
robust  or  alternatively,  for  given  control  gain,  can  be  used  to  determine  the 
range  of  the  size  of  allowable  perturbations  for  stability.  A  VTOL  aircraft 
example  [23]  is  considered  in  which  variations  of  parameters  due  to  changes 
in  air  speed  are  taken  as  perturbations  and  a  constant  gain  linear  state 
feedback  control  law  is  presented  to  attain  stability  robustness  for  these 
perturbations.  Various  design  implications  are  discussed  based  on  this  example. 
3.1  Linear  State  Feedback  Control  Design  Using  Perturbation  Bound  Analysis 
As  before,  consider  the  linear,  time  invariant  system  described  by 

£  =  Ax  +  Bu  (3 . 1) 

y  =  Cx 

where  x  is  nxl  state  vector,  the  control  u  is  mxl  and  output  y  (the  variables 
we  wish  to  control)  is  kxl.  The  matrix  triple  (A,B,C)  is  assumed  to  be 
completely  controllable  and  observable.  Let  the  control  law  be  given  by 
u  =  Gx. 
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For  this  case,  the  nominal  closed  loop  system  matrix  is  given  by 

A  =  A  +  BG  £3  < 

and  G  is  such  that  A  is  asymptotically  stable. 

Let  AA,  AB  and  AG  be  the  perturbation  matrices  formed  by  the  maximum 
modulus  deviations  expected  in  the  individual  elements  of  matrices  A,  B 
and  G  respectively.  Then  for  this  structured  perturbation  case,  the 
perturbed  system  matrix  of  (3.1)  is  given  by 


A^  =  (A  +  AA)  +  (B  +  AB)  (G  +  AG) 


(3.3) 


Design  Observation  3.1:  The  perturbed  system  matrix  of  (3.3)  i 
stable  for  all  perturbations  bounded  by  AA,  AB  and  AG,  if 


is 


£==£..  =  [AA  +  AB  G  +  (B  +  AB)  AG]  .  .  <  - 7= - —  - 

ij  max  m  m  lj  max  a  (P  U  )  =  U 

max  m  e  s  Y 


(3.4) 


where  P  satisfies  PA  +  A  P  +  21  =0  and  P  =  Ip  I. 

n  mij  1  ij  1 

Note  that  the  'error'  matrix  in  the  left  hand  side  (l.h.s.)  of  (3.4) 
comprises  the  perturbation  matrices  AA,  AB  and  AG  as  well  as  nominal  matrices 
G  and  B.  In  practice,  the  following  situations  may  arise. 

Case  A:  The  matrices  AA,  AB  and  AG  are  known  (of  course  G  and  B  are  known 

m 

too):  for  this  case,  the  matrix  in  the  right  hand  side  (r.h.s.)  is  simply 

formed  by 


U  .  .  =  e .  . /z  ( q  cr'v 

ei-j  ij 

Case  B:  Matrices  AA,  AB  and  AG  are  partially  known:  in  this  case  one  needs 


at  least  any  two  ratios  among  (AA. .  /AB. .  )  or  (AA  /AG  ) 

13  max  13  max  ij  max  ij  max 

or  (AB-  •  „__/AG.  .  )  to  be  able  to  form  the  U  matrix  (in  which  case  the 

13  max  jlj  max  0 

Ue  matrix  contains  some  unity  elements  at  the  appropriate  places). 

Case  C;  No  information  on  the  perturbation  matrices  in  the  l.h.s.  is  known: 


for  this  case,  one  has  to  go  for  the  worst  case  situation  and  replace  U  by 


U  •  In  other  words  the  stability  condition  becomes 
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e  <  a  (P  u  ) 
max  m  n  s 


(3.6a) 


where  E  =  [AA  +  AB  G^  +  (B  +  AB)  AG].  (3.6b 

It  may  be  seen  that  similar  development  applies  to  the  Unstructured1 
perturbation  case  also  and  hence  the  details  are  not  given  here. 

3.1.1  Stability  Robustness  Index  and  Design  Algorithm 

A  first  glance  at  the  stability  condition  (3.4),  may  suggest  that  the 
bound  (for  structured  perturbations)  can  serve  as  a  measure  of  stability 
robustness  and  that  higher  this  bound,  better  is  the  stability  robustness. 
It  may  thus  be  concluded  that  the  control  gain  G  which  makes  y^  maximum  is 
the  gain  to  look  for.  However  a  closer  look  at  the  condition  reveals  that 
both  l.h.s.  and  r.h.s.  terms  of  the  conditions  (3.4)  are  functions  of  the 
control  gain  G  and  a  variety  of  control  gains  G  may  satisfy  the  proposed 
stability  condition  for  given  perturbations.  It  is  thus  more  appropriate 
to  consider  the  relative  difference  of  these  two  terms  in  comparing  and 
synthesizing  different  controllers  from  stability  robustness  point  of  view. 
To  this  end  we  define,  as  a  measure  of  stability  robustness,  an  index 
called  'Stability  Robustness  Index*  given  by 


3  =  (y  -  e)  (for  structured  perturbations) 

S  •  R.  Y 


(3.7a) 


(3.7b) 


for  Case  A  situation  and 

B  =  y  (structured  perturbation) 

S  •  R.  Y 

for  Cases  B  and  C  situations  described  before. 

By  this  definition,  for  Case  A,  8  >0  corresponds  to  the  stability 

robustness  region. 

Even  though  the  'stability  robustness  index*  defined  above  plays  the  role 
of  a  'stability  margin',  it  is  quite  different  in  its  interpretation  from  the 
standard  notions  of  'stability  margin*.  One  classical  relative  stability 
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margin  in  Single  Input  Single  Output  (SISO)  stable  systems  is  the  distance, 
a,  of  the  nearest  characteristic  root  from  the  imaginary  axis.  Another  margin 
of  stability  in  SISO  systems  is  the  standard  gain/phase  margin ,  namely  the 
gain/phase  changes  that  can  be  tolerated  before  the  system  becomes  unstable. 
There  are,  of  course,  investigations  being  carried  out  to  generalize  these 
notions  for  multivariable  systems  [1].  In  contrast  to  these  interpretations, 
6  is  a  stability  margin  given  in  terms  of  the  tolerable  system  parameter 

perturbations.  That  is,  higher  the  index  3  ,  higher  is  the  tolerable 

S .  R. 

parameter  perturbation  for  stability.  Clearly  these  different  notions  of 
stability  margins  do  not  imply  one  another,  except  in  special  circumstances. 
Thus  higher  $  does  not  necessarily  mean  higher  a  and  vice  versa.  It  may 

also  be  noted  from  (3.4)  that  for  linear  systems  when  AA  =  0,  AB  =  0  and  AG  ^  0 
(unknown),  gg  R  (=py )  may  be  regarded  as  a  'gain  margin',  albeit  in  a 
restricted  sense. 


3 >1.2  Control  Design  Algorithm  for  Robust  Stability 

It  is  clear  from  the  expressions  for  (3.4),  the  'error  matrix'  (3.6b) 
and  3  (3.7)  that  these  quantities  depend  on  the  control  gain  G  and  as  gain 

G  is  varied  g  changes.  In  order  to  plot  the  relationship  between  g 

S.R.  S.R. 

and  the  gain  G,  we  need  a  scalar  quantitative  measure  of  G.  For  this  we  can 
either  use 


J 

cn 


s 


a 

max 


(G) 


(3.8a) 


or 


CO  oo 

J  =  [/  (uTu)dt]1/2  =  [/  xTGTGxdt] 1/2 

cn  J o  J o 

where  J  denotes  a  measure  of  'nominal  control  effort*, 
cn 

The  variation  of  3  with  the  control  effort  J 

S.R.  cn 

on  the  perturbation  matrices  and  on  the  behavior  of  the 


(3.8b) 

is  very  much  dependent 
Lyapunov  solution. 
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which  cannot  be  described  analytically  in  a  straightforward  way.  Assuming 
stability  robustness  is  the  only  design  objective,  the  design  algorithm 
basically  consists  of  picking  a  control  gain  that  maximizes  stability  robust¬ 
ness  (Sg  R  ) •  Specifically  the  algorithm  involves  determining  the  index 
^S.R.  and  the  contro1  effort  J  for  different  values  of  the  control  gain  G 
and  plotting  these  curves.  These  design  curves  can  then  be  used  to  pick  a 
gain  that  achieves  a  high  R  •  The  algorithm  thus  provides  a  simple  constant 
gain  state  feedback  control  law  that  is  robust  from  stability  point  of  view. 

The  algorithm,  for  given  perturbations,  can  be  used  for  selecting  the  range 
of  control  effort  for  which  the  system  is  stability  robust  or  alternatively 
for  given  control  effort,  can  be  used  to  determine  the  range  of  allowable 
perturbations  for  stability. 

In  the  next  section  a  widely  used  VTOL  aircraft  control  problem  with 
varying  flight  conditions  is  considered  and  a  stability  robust  controller 
is  synthesized  using  the  proposed  methodology. 

3.1.3  Robust  Control  Design  for  VTOL  Aircraft 

In  this  section,  we  recall  the  VTOL  aircraft  (helicopter)  problem 
described  in  section  (2.3). 

The  linearized  model  of  the  VTOL  aircraft  in  the  vertical  plane  is 
described  by 

x  =  (A  +  AA) x  +  (B  +  AB)u,  x (0)  =  XQ  (3.9) 

The  components  of  the  state  vector  xeR4  and  the  control  vector  ueR2  are 
given  by 

->  horizontal  velocity  (knots) 
x2  ->  vertical  velocity  (knots) 
pitch  rate  (degree/s) 
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x4  pitch  angle  (degrees) 

->  'collective'  pitch  control 
u2  ^  1  longitudinal  cyclic'  pitch  control 

Essentially,  control  is  achieved  by  varying  the  angle  of  attack  with 
respect  to  air  of  the  rotor  blades.  The  collective  control  U;L  is  mainly 
used  for  controlling  the  motion  of  the  aircraft  vertically  up  and  down. 
Control  u^  is  basically  used  to  control  the  horizontal  velocity  of  the 
helicopter. 


In  [23] ,  the  mathematical  model  is  presented  assuming  the  nominal 

airspeed  to  be  135  knots.  For  this  nominal  case  (A  A  and  AB  =  0)  the  matrices 
A  and  B  are  given  by 


A  = 


-0.0366 

0.0482 

0.1002 

0 


0.0271 

-1.01 

0.3681 

0 


0.0188 

0.0024 

-0.707 

1 


-0.4555 

-4.0208 

1.42 

0 


B  = 


0.4422  0.1761 

3.5446  -7.5922 

-5.52  4.49 


_0  0 

We  take  the  initial  condition  to  be 
T 

x  (0)  =  [0.85  0.15  0  -0.05] 


(3.10a) 


(3.10b) 


(3.11) 


As  the  airspeed  changes  significant  changes  take  place  in  the  elements 

a32'  a34  and  b21*  let  us  consider/  for  illustration  purposes,  the  following 
range  of  variations  in  these  parameters. 

Case_l:  0.3545  £  a 32  =  0.3681  <_  0.3817 

1-31  <_  a34  =  1.42  _<  1.53 
3.39  <_  b  =  3.544  <_  3.702 
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(3.12) 


i.e.,  lAa32l  =  °*0136'*  lAa34l  =  °-1:L'  I  Ab2i  I  =  °-157 

Case  2:  |Aa32|  =  0.041;  |Aa34|  =  0.332;  |m>21|  =  0.47  (3.13) 

Case  3:  lAa32l  =  0.068;  |  Aa34 1  =  0.553;  |  At>  |  =  0.78  (3.14) 

Case  4:  1 Aa32 1  =  0.095;  j  Aa34 1  =  0.774;  |Ab21|  =  1.097  (3.15) 

Case  5;  |  A  a3  2 1  =  0.1363;  |Aa3J  =  1.106;  |  At>21|  =  1.5674  (3.16) 


In  other  words,  the  matrices  AA  and  Ab  are  known  for  these  five  cases. 
For  designing  the  nominal  state  feedback  control  that  stabilizes  the 
nominal  closed  loop  system,  for  this  example,  we  employ  the  standard  linear 
quadratic  optimal  control  design  algorithm  assuming  a  performance  index 


CO 

r  T  T 

J  =  Jo  (x  Qx  +  p^u  R^u)dt 


where  Q  and  Rq  are  (nxn) ,  and  (mxm)  symmetric  positive  definite  matrices  and 
Pc  is  a  scalar  variable  used  for  designing  the  control  gain  G. 

For  this  case,  the  nominal  closed  loop  system  matrix  is  given  by 


A  =  A  +  BG, 


BTK/Pc 


(3.17) 


and 


T 

KA  +  A  K  -  KB 


T  T 

•BK+CQC=0 


and  A  is  asymptotically  stable.  For  each  case  using 


(3.18) 


Q  = 


0,04 

0.01 

0 

0.01 

0.01 

0.025 

0.01 

0 

and  R  =  I„  and  p 

0 

0.01 

0.01 

0 

o  2  c 

0.01 

0 

0 

0.01 

design 


variable,  the  standard  optimal  IQ  regulator  control  gain  and  the  corresponding 

control  effort  J  of  (3 . 8b)  are  computed  and  8  is  calculated.  The  plots  of 

cn  S  •  R. 

3  „  vs.  J  for  these  five  cases  are  shown  in  Fig.  4. 

►d  .  cn 


-54- 


In  interpreting  these  plots,  it  is  to  be  recalled  that  the  region  of 

control  effort  for  stability  robustness  is  the  region  in  which  3  >  0. 

S  •  R. 

Design  observation  3 .1':  From  these  plots,  it  may  be  observed  that  as  the 
parameter  perturbation  range  is  increased,  the  range  of  control  effort  for 
stability  robustness  decreases. 

Design  observation  3-. 2:  For  given  set  of  parameter  perturbations,  there  is  a 
unique  control  effort  for  which  B  is  maximum.  Obviously  this  is  the 

control  effort  (and  the  corresponding  control  gain)  that  we  seek,  assuming 
practical  constraints  are  satisfied. 

Design  observation  3.3:  It  is  to  be  noted  that  for  all  these  cases,  the 

maximum  0  occurs  almost  at  the  same  control  effort.  However  it  is  also 

b  •  K. 

to  be  kept  in  mind  that  the  range  of  variations  considered  in  cases  2  through 
5  are  simply  some  multiples  of  the  range  of  case  1. 


Now  let  us  consider  the  cases  when  there  are  parameter  perturbation  in 
only  one  of  the  matrices,  A  or  B.  Accordingly,  we  consider  the  following 
ranges . 


Case 

6: 

Aa32 

|  =  0.3018; 

1 ia34 1 

=  1.300; 

KJ  - 0 

Case 

7: 

lAa32 

|  =  0.1366; 

l4a34i 

=  1.106; 

KJ  - 0 

Case 

8: 

lAa32l 

II 

o 

SI 

o 

CD 

OJ 

1  - 

|Ab21l  - 

2.5671 

Case 

9: 

lAa32l 

-  0;  |Aa34 

|  -  0; 

1.5674 

Figure  5  corresponds  to  these  cases.  It  may  be  observed  from  these 
plots  that,  as  before,  it  turns  out  that  smaller  the  size  of  the  perturbation, 
the  more  is  the  control  effort  range.  But  consider  cases  6  and  9.  From 
plots  corresponding  to  these  two  cases,  it  can  be  seen  that  the  control  range 
(for  stability)  for  perturbations  in  A  is  larcrer  than  the  control 
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range  for  perturbations  in  B  indicating  variations  in  matrix  B  are 
more  critical  from  stability  robustness  point  of  view. 

Extention  to  LQG  Regulators  and  Dynamic  Compensator  Design 

Clearly  the  proposed  methodology  can  be  extended  to  Linear  Quadratic 
Gaussian  (LQG)  regulators  and  dynamic  compensators  in  a  straightforward  way. 
However,  the  design  implications  in  this  case  could  be  different.  This 
extension  is  suggested  as  a  future  research  topic. 


J  (AA  /  0,  AB  ^  0). 


J  (AA  =  0,  or  AB  =  0) 
cn 
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IV.  ANALYSIS  FOR  PERFORMANCE  ROBUSTNESS  OF  LINEAR  REGULATORS 


At  this  stage,  it  is  instructive  to  pause  for  a  moment  and  recapitulate 
the  material  that  has  been  discussed  in  the  previous  sections.  The  basic 
aspect  addressed  was  that  of  'stability  robustness'  namely  the  assumption 
that  the  fundamental  'acceptable1  behavior  is  the  stability  of  the  closed 
loop  system  in  the  presence  of  perturbations.  However,  in  many  practical 
situations,  guaranteeing  stability  alone  may  not  be  sufficient  to  carry  out 
the  mission  objective.  This  is  particularly  so  in  the  'linear  regulator* 
problems  with  stringent  requirements  on  the  'regulation*  aspect,  especially 
in  aerospace  applications. 

It  is  the  purpose  of  the  control  design  to  'regulate'  the  controlled 
variables  (or  state)  and  the  performance  of  the  regulator  is  deemed 
'acceptable'  if  the  'closed  loop  regulation  cost'  in  the  presence  of  these 
perturbations  is  within  a  bound  specified  by  the  designer  or  the  design 
specifications.  Thus,  not  only  stability  but  a  specified  degree  of  per¬ 
formance  (regulation)  are  required.  This  aspect  is  termed  the  'Performance 
(Regulation)  Robustness'  and  this  section  addresses  the  aspect  of  analysis 
and  design  of  linear  quadratic  regulators  for  this  kind  of  'performance 
robustness 1 . 

4. 1  Brief  Review  of  Literature 

The  ' Robustness '  of  linear  multivariable  regulators  has  been  analyzed 
by  many  researchers.  Anderson  and  Moore [38]  have  shown  that  the  single  input 
Linear  Quadratic  State  Feedback  (IQ)SF)  designs  have  a  phase  margin  of  60° 
and  infinite  gain  margin.  Recent  studies  by  Wong  and  Athans  [39]  and  Safonov 
and  Athans  [40]  conclude  that  the  multi  input  LQ  regulators  possess  a 
guaranteed  gain  margin  -  6  db  to  +<*>db  and  a  phase  margin  of  +60°  in  all 
channels.  It  is  shown  by  them  that  no  such  guaranteed  margins  can  be  given 
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for  Linear  Quadratic  Gaussian  (LQG)  regulators.  The  robustness  results 
presented  in  these  above  references  are  primarily  concerned  with  allowable 
perturbation  in  the  feedback  gains.  However  these  robustness  properties 
cannot  directly  be  extended  to  the  case  of  perturbations  in  system  matrices 
(i.e.,  parameter  variations).  In  fact,  it  is  shown  recently  by  Soroka  and 
Shaked  [41]  that  despite  possessing  attractive  gain  and  phase  margin 
properties,  the  LQ  regulator  suffers  from  poor  stability  robustness  from 
parameter  variations  point  of  view.  It  is  for  this  reason  that  we  focus  our 
attention  on  the  analysis  and  design  for  robustness  with  respect  to  parameter 
variations . 

Ly  and  Cannon  [42]  give  stochastic  representation  for  parameter  uncertainty 
whereas  this  report  considers  the  uncertain  parameters  to  vary  within  certain 
bounds.  The  work  of  Leitmann  [24],  Barmish  and  colleagues  [26]  on  ultimate 
boundedness  control  of  mismatched  uncertain  system  addresses  the  core 
problem  of  regulation  but  the  control  law  given  is  nonlinear  while  the  aim 
of  this  section  is  to  utilize  the  simplicity  of  linear  feedback  control  law  of 
the  standard  optimal  LQ  regulators.  Lunze  [6] ,  using  comparison  systems 
concept  and  Owens  and  Chotai  [7]  by  combining  frequency  domain  and  time  domain 
treatments,  present  attractive  schemes  involving  more  general  perturbations. 
Their  methods  do  not  specifically  address  regulators  with  quadratic  performance 
indices.  Desoer  et  al  [8]  have  established  conditions  for  stability  robust¬ 
ness  of  linear  multivariable  interconnected  systems  for  sufficiently  small 
perturbations . 

But  these  methods  do  not  specifically  address  the  type  of  ‘regulation 
robustness'  outlined  above.  In  the  context  of  regulation  robustness,  Fujii 
and  Mizushima  [43]  address  the  robustness  of  the  optimality  property  of 
linear  regulators  but  the  analysis  is  restricted  to  small  parameter  variations 
and  is  carried  out  in  frequency  domain.  Rolnik  and  Horowitz  [44]  present  a 


-58- 


method  suitable  for  large  parameter  variations,  again  in  frequency  domain, 
but  for  a  third  order  system.  In  time  domain,  interesting  work  is  reported  on 
the  robustness  of  linear  regulators  by  Sezer  and  Siljak  [45,46,47]  using  the 
concept  of  suboptimality  index  which  is  used  as  a  measure  for  system  robust¬ 
ness  both  with  respect  to  parameter  variations  as  well  as  suboptimal  control 
law.  The  analysis  is  done  for  time  invariant  perturbations  in  the  context  of 
interconnected  systems.  Krishnan  and  Brzezowski  [48]  present  an  iterative 
design  algorithm  such  that  prescribed  trajectory  insensitivity  is  achieved 
but  the  method  is  aimed  at  time  invariant  perturbations.  Barnett  and  Story 
[49]  consider  a  specific  variation  for  which  the  perturbed  cost  is  same  as 
the  nominal  optimal  cost.  McClamroch  et  al  [50]  offer  a  method  suitable 
only  for  small  number  of  uncertain  parameters  and  fow  low  order  systems.  An 
interesting  method  by  Burghardt  [51]  proposes  the  weighting  matrices  to  be 
functions  of  the  uncertain  parameters.  The  method  is  ad-hoc  and  iterative 
in  nature.  The  method  of  Evans  and  Xyaoni  [52]  considers  time  invariant 
perturbations  and  employs  parameter  optimization  techniques  to  design  the 
controller. 

From  a  different  perspective,  there  is  interesting  literature  on  the 
analysis  and  design  procedures  for  tolerable  perturbations  (Perturbation 
Bound  Analysis)  for  robust  stability  and  regulation.  References  [9] -[15] 
deal  with  bounds  for  robust  stability  alone  with  no  concern  for  the  regulation 
robustness.  Rissanen  [53]  presents  bounds  on  the  linear  perturbation  of  the 
nominal  regulator  such  that  the  perturbed  regulation  cost  is  below  the  nominal 
cost  times  a  given  constant  greater  than  1.  The  proposed  bound  holds  good 
for  time  varying  perturbations  since  the  method  uses  Lyapunov  approach  but 
no  synthesis  procedure  is  given.  Sarma  and  Deekshatulu  [54]  address  a  similar 
problem  but  the  method  is  suitable  only  for  small  number  of  uncertain  para¬ 
meters  and  low  order  systems.  Perhaps,  one  of  the  best  known  design  methods 
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in  time  domain  for  time  varying  perturbations  is  the  method  of  Guaranteed 
Cost  Control  of  Chang  and  Peng  [13}.  (Ref.  [55]  by  Vinkler  and  Wood  extends 
this  method  for  time  invariant  perturbations.)  The  design  procedure  proposed 
in  this  paper  belongs,  in  principle,  to  the  framework  of  Guaranteed  Cost 
Control  method  m  the  sense  that  it  is  a  design  method  in  time  domain  suitable 
for  time  varying  perturbations,  belonging  to  bounded  intervals.  The  relative 
merits  of  these  methods  are  discussed  later.  The  proposed  procedure  basically 
formulates  the  available  Perturbation  Bound  Analysis  into  an  algorithm  to 
design  controllers  for  regulation  robustness.  Towards  this  direction,  in  this 
section,  we  present  perturbation  bounds  for  'robust  regulation'  and  then  in 
Section  V  present  a  design  algorithm  to  synthesize  controllers  for  robust 


regulation . 

4-2  Perturbation  Bounds  for  Robust  Regulation 

In  this  section,  we  briefly  review  the  upper  bounds  for  robust  stability, 
presented  in  previous  sections  (for  structured  uncertainty) . 

Upper  Bounds  for  Robust  Stability:  Consider  x  =  Ax  (41) 

X  =  [A  +  E  (t)  ]  X 


The  system  of  (4.1)  is  stable  if 


a  .  (Q) 

min  * 

6  < -  ■  ■■  =  n 

0  (P  U  )s  y 

max  me  1 


(4.3a) 


(4.3b) 


where  P  satisfies  the  equation 

PA  +  ATP  +  2Q  =  0 

and  Ug  is  an  nxn  matrix  whose  entries  are  such  that 

U  .  .  =  e  /e  (0  <  u  .  .  <  1) 
eil  13  —  ei]  — 

Note  that  U^  (and  thus  u  )  can  be  computed  even  if  one  knows  only  the  ratio 

instead  of  knowing  e. .  (and  e)  separately. 

J  1 1 
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One  important  difference  in  the  formula  for  y^  of  (3.4)  above  and  the 
bound  presented  in  (4.3)  is  the  presence  of  am£n(Q)  in  the  numerator  of  (4.3). 
However,  u  of  (4.3)  is  the  more  general  expression  (Refs.  [13] -[53])  and  u 

y  y 

of  (3.4)  is  simply  a  special  case  of  (4.3)  where  Q  =  I.  The  reason  we  now 
use  y^  of  (4.3)  is  the  fact  that  in  this  section  we  intend  to  relate  this 
bound  for  robust  stability  to  the  bound  we  are  going  to  present  for  robust 
regulation.  When  only  robust  stability  is  the  primary  concern  (with  no  concern 
for  robust  regulation)  then,  of  course,  we  can  use  y  of  (3.4)  since  0=1 

y 

maximizes  this  bound  as  shown  in  [14] . 

We  now  define  precisely  our  interpretation  of  1  robust  regulation*  and 
then  present  perturbation  bounds  to  achieve  this  kind  of  ’robust  regulation1 
Upper  Bounds  for  Robust  Regulation 
Let 


a  oo  rr 

J  =  /  X  Q  x  dt 

n  o 


(4.4) 


A  T 

J  =  /  X  Q  Xdt 
P  o 

We  deem  the  performance  of  the  above  system  to  be  'acceptable'  if 

J  <  J 
p  m 

where  is  the  ’maximum  tolerable  regulation  cost'.  Suppose 


(4.5) 


(4.6) 


J  =  k  J 
m  n 

P.e^Pit:*-on :  The  system  of  (4.2)  is  ’robust*  from  performance  (regulation) 

point  of  view  if 


J 

P 


<  k  J  (=  J  ) 
n  m 


for  some  given  k  >  1. 


(4.8) 


With  this  definition ,  the  bounds  presented  by  Rissanen  [53]  on  E  can 
be  extended  as  follows. 
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i)  For  unstructured  perturbations:  For  system  of  (4.2),  the  performance 
robustness  criterion  of  (4.8)  is  satisfied  if 


v-i  ^^(Q) 

a  (E  (t) )  <  -2i2 - 

max  —  k 


a  (P) 
max 


(4.9) 


where  P  satisfies  the  equation  (4.3b). 

ii)  For  structured  perturbation:  For  system  of  (4.2),  the  performance 
robustness  criterion  of  (4.8)  is  satisfied  if 


e  < 


k-1  °min(Q) 


—  k  a  (P  U  ) 
max  m  e  s 


(4.10a) 

(4.10b) 


Henceforth,  the  discussion  will  be  focused  on  'structured*  perturbation 
case,  keeping  in  mind  that  similar  development  applies  to  the  'unstructured' 
perturbation  case  also. 

From  equations  (4.3),  (4.10)  it  is  clear  that  the  perturbation  bound  for 
'robust  regulation'  is  always  less  than  or  equal  to  the  bound  for  'robust 
stability ' . 

Extension  to  LQ  Regulators 

Consider  the  linear,  time  invariant  system  described  by 

x  =  Ax  +  Bu  (4.11) 

y  =  Cx  (C  is  full  rank)  (4.12) 

where  x  is  nxl  state  vector,  the  control  us  is  mxl  and  output  y  (the  variables 
we  wish  to  control)  is  kxl .  Let  the  performance  index  for  the  above  system  be 
given  by 

CJO  *p 

J  =  fQ  (y  Qy  +  PQU  RQu)dt  (4.13) 

where  Q  and  are  (kxk)  and  (mxm)  symmetric  positive  definite  matrices, 
respectively,  and  is  a  scalar  variable  used  for  design.  The  matrix  triple 
(A , B , C)  is  assumed  to  be  completely  controllable  and  observable. 
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For  this  case ,  the  nominal  closed  loop  system  matrix  is  given  by 

*  =  A*  (4.14a) 

(4.14b) 


A  =  A  +  BG,  G  =  -R  1  BTK/ 

o  c 


and 


R 


-1 


T  o  T  T 

KA  +  A  K  -  KB  -  B  K  +  C  QC  =  0 


(4.14c) 


and  A  is  asymptotically  stable. 

Let  AA,  AB  and  AG  be  the  perturbation  matrices  formed  by  the  maximum 
modulus  deviations  expected  in  the  individual  elements  of  matrices  A,  B  and 
G,  respectively.  Then  for  this  structured  perturbation  case,  the  perturbed 
system  matrix  is  given  by 

X  =  X  (4.15a) 

where 

Ap  =  (A  +  AA)  +  (B  +  AB)  (G  +  AG)  (4.15b) 

The  basic  premise  of  Regulation  Robustness1  involves  the  definition  of 
acceptable  performance .  In  linear  regulator  problems ,  it  is  the  purpose  of 
the  control  design  to  'regulate'  the  state  and  the  efficacy  of  the  design 
and  is  measured  by  the  'closed  loop  regulation  cost'  defined  by 

(4.16) 

subject  to  equation  (4.1).  But  in  the  presence  of  perturbations,  with  the 


A  oo  —  _  m 

Jxn  =  fQ  x  Qx  dt.  Q  =  c  QC 


fixed  controller  of  (4.14),  the  'perturbed  regulation  cost'  defined  by 

A  .«>  T — 

J  =  /  X  QX  dt 
xp  o 

subject  to  equation  (4.15),  may  deviate  considerably  from  J  of  (4  16) 

xn 

However,  the  design  specifications  may  dictate  that  the  performance  is 
’unacceptable'  if  the  'closed  loop  regulation  cost'  exceeds  a  specific  given 
value.  Thus,  we  deem  the  performance  of  the  regulator  to  be  'acceptable'  if 
JXp  <  Jm  course  Jxn  <  Jm)  where  is  the  'maximum  tolerable  regulation 
cost'  (assuming  the  control  effort  to  achieve  this  is  acceptable) . 
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As  before,  letting 


J  -  k  J  ,  k  >  1 
m  xn 

(4.18) 

we  say,  the  linear  quadratic  regulator  is  Robust'  from  performance 


(regulation)  point  of  view  if 

J  <  k  J 

xp  xn 

(4.19) 

for  some  k  >  1 . 

Design  observation  4,1a;  The  perturbed  L Q  Regulator  system  matrix  of  (4.15) 


is  stable  for  all  perturbations  bounded  by  A,  B,  and  G  if 

E  =  [AA  +  AB  G  +  (B  +  AB)  G]  .  . 

m  m  13  max 

a  .  (CTQC) 
mm 

<  - —  ~  ij 

a  (P  u  )  y 

max  m  e  s 

(4.20a) 

—  — t  T 

where  P  satisfies  PA  +  A  P  +  2C  QC  =  0 

(4.20b) 

Design  observation  4.1b:  The  perturbed  L Q  Regulator  satisfies  the  regulation 

robustness  criterion  of  J  <  k  J  =  J  if 

xp  xn  m 


k-1 

e  <  y  -  y 

k  y  yp 

(4.21) 

where  e  and  y^  are  as  defined  in  (4.20). 

Having  developed  the  perturbation  bounds,  we  now  present  a  design 
algorithm,  using  these  bounds,  to  synthesize  controllers  for  robust 
regulation . 
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V.  SYNTHESIS  OF  CONTROLLERS  FOR  ROBUST  REGULATION 


In  Section  III,  a  synthesis  procedure  for  linear  systems  to  have  stability 
robustness  was  presented  using  the  perturbation  bounds  for  robust  stability.  In 
that  discussion,  it  was  assumed  that  the  main  design  objective  was  stability 
alone  with  no  concern  for  regulation.  In  this  section,  we  use  essentially  the 
same  principle  to  synthesize  controllers  for  linear  regulators  for  maintaining 
not  only  stability  but  an  'acceptable'  level  of  performance  (or  regulation). 
Even  though  the  concept  is  similar,  the  design  implications  are  somewhat 
different  in  this  proposed  algorithm. 

Introducing  a  quantitative  measure  called  'Performance  Robustness  Index'  a 
design  algorithm  is  presented  by  which  one  can  achieve  a  trade  off  between 
nominal  performance  and  regulation  robustness.  The  proposed  method  is 
illustrated  with  the  help  of  a  simple  example  and  results  are  discussed. 

5.1  Performance  Robustness  Index  and  Design  Algorithm 

From  (4.18)  and  (4.21),  it  can  be  seen  that  the  performance  robustness 
criterion  involves  two  quantities,  namely  Jm,  the  tolerable  regulation  cost  and 
Uyp'  the  tolerable  perturbation.  The  tolerable  regulation  cost  Jm  can  be 
looked  at  from  two  viewpoints:  either  i)  it  can  be  specified  as  a  constant 
number  in  which  case  the  index  k  and  the  nominal  regulation  cost  Jxn  are  such 
that  their  product  yields  the  specified  constant  number,  or  ii )  it  can  be  given 
in  terms  of  a  constant  k  times  the  nominal  regulation  cost  Jxn.  In  the  former 
case  Jm  is  constant  whereas  in  the  latter  case  Jm  is  an  explicit  function  of 
Jxn*  In  either  case  it  may  be  noted  that  both  the  left  hand  side  (l.h.s.) 
terms  of  the  condition  are  functions  of  the  control  gain  G  and  a  variety  of 
control  gains  may  satisfy  the  proposed  performance  robustness  condition  for 
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given  perturbations  and  Jm.  Hence,  to  compare  and  synthesize  different 
controllers  from  robustness  point  of  view,  there  is  a  need  for  defining  a 
’measure'  of  performance  (regulation)  robustness.  From  (4.18  ),  (4.21)  it  is 
clear  that  regulation  robustness  may  be  measured  either  in  terms  of  the 
perturbation  bound  ^yp)  or  regulation  bound  (Jm),  depending  upon  the 
specifications.  Accordingly,  we  define  a  generic  index  called  'Performance 
Robustness  Index,  1  as  a  measure  of  performance  (regulation)  robustness. 

In  what  follows,  we  delineate  the  different  situations  one  encounters  in  a 
design  environment  and  the  corresponding  interpretation  for  8p.R# .  For  each 
case  a  design  algorithm  is  presented  based  on  this  index. 

Case  ia:  (constant)  is  given  (i.e.  k  is  known)  and  e  is  known  a  priori 

(i.e.  AA,  AB  ...  are  known). 

Jm  (variable)  is  known  (i.e.  constant  k  is  given)  and  e  is  known  a 
priori  (i.e.  AA,  AB  ...  are  known). 

Since  in  these  two  cases,  the  tolerable  regulation  Jm  is  specified,  we 
take  the  tolerable  perturbation  ^yp  as  the  design  objective.  Accordingly,  one 
appropriate  measure  of  performance  robustness  is  the  relative  size  of  the 
residual  perturbation.  Thus  we  define 

Bp.R.  =  (Myp  -e )  (5.1 ) 

Since  s  is  presumed  known,  one  has  to  check  whether  the  given  perturbation 
satisfies  the  condition  (4.21 )  or  not.  As  e  and  ^yp  are  both  functions  of  the 
control  gain  G,  the  performance  robustness  region  is  given  by  8P  R  >  0. 

Design  Algorithm;  For  this  case,  the  design  algorithm  is  as  follows: 

1 )  ^-specific  requirements  on  the  nominal  performance  (i.e.  on  Jxn) : 

For  this  situation,  the  algorithm  basically  consists  of  picking  a  control 
gain  (i.e.  pc)  that  maximizes  Sp.R.  (=  ^yp“e  )  in  ^rie  region  8p.R.  >  0. 
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2 )  Concern  for  both  nominal  performance  as_  well  as  performance  robustness; 

For  this  situation,  the  algorithm  consists  of  plotting  3p  R  vs.  the 

A  00 

nominal  control  effort  Jen  |  (=  /  uTu  dt)1/2  and  Jxn  vs.  Jcn  and  picking  a 

o 

control  gain  that  achieves  a  reasonable  trade  off  between  these  variables. 

Note  that  these  design  curves  can  be  used,  for  given  set  of  perturbations, 
to  select  the  range  of  nominal  control  effort  for  which  the  system  is 
performance  robust. 

Case  ii  JLsJL*  Jm  <c  onstant)  is  given  (i.e.  k  is  known)  and  e  is  not  known. 

Case  ii(b) :  (variable)  is  known  (i.e.  constant  k  is  given)  and  e  is  not 

known . 

This  case  is  essentially  the  same  as  before  except  that,  this  time  we 
define 


3 


P.R . 


y 


yp 


(5.2) 


In  other  words,  tolerable  perturbation  is  still  taken  as  the  design  objective. 
Design  Algorithm:  Essentially  the  same  as  before  except  that  we  do  not  need 
the  requirement  of  checking  whether  $p#R#  is  positive  or  not  (it  is  always 
positive ) . 

Note  that  for  this  cases,  the  design  curves  can  be  used,  for  given  control 
effort,  to  determine  the  range  of  size  of  allowable  perturbation. 

For  case  (i)  and  (ii),  the  performance  robustness  index  is  in  terms 

of  tolerable  perturbation  bound  and  thus  higher  indicates  better 

performance  robustness. 

Case  iii ( a) :  is  now  known  (i.e.  k  is  not  given  and  e  is  known). 

Since,  in  this  case,  the  tolerable  regulation  cost  Jm  (and  thus  k)  is  not 
given,  we  will  use  it  as  a  design  objective.  In  other  words  our  aim  would  be 
to  minimize  Jm  under  the  worst  case  situation  of  the  given  perturbation  being 
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equal  to  the  bound  Pyp.  Accordingly,  from  the  condition  (4.21),  we  write 


(5.3) 


from  which  we  get 


y 


k 


(5.4) 


We  now  define  the  performance  robustness  index  as 


(5.5) 


where  k  is  given  by  (5.4)  and  Jm  is  required  to  be  positive.  With  this 
definition,  is  expressed  in  terms  of  the  tolerable  regulation  bound 


(rather  than  the  perturbation  bound  as  in  cases  i  and  ii).  The  above 
definition  is  such  that  lower  the  value  of  Jm,  better  (or  higher)  is  the 
performance  robustness. 

Evidently,  for  this  situation  one  would  like  to  achieve  a  low  value  of  Jm 
( >0 )  from  high  performnce  robustness  point  of  view  and  a  low  value  of  Jxn  from 
nominal  performance  point  of  view  (assuming  adequate  control  effort  is 
available).  It  is  important  to  note  that  Hy,  e  and  Jxn  are  all  functions  of 
the  control  gain  G  (or  the  control  weighting  factor  Pc )  and  their  variation 
with  G  (or  Pc)  is  very  much  dependent  on  the  perturbation  matrices  and  on  the 
behavior  of  the  Lyapunov  solution  which  cannot  be  ascertained  analytically  or 
qualitatively  in  a  straightforward  manner,  although  it  may  be  possible  to  make 
general  comments  in  some  special  cases  of  perturbation.  Hence  in  some  problems 
it  may  be  possible  to  achieve  both  high  (low  Jm)  as  well  as  low  Jxn  with 

the  same  control  gain  whereas  in  some  cases,  a  trade  off  between  St,  D  and 
may  be  necessary. 

Design  Algorithm:  For  this  case,  the  design  algorithm  consists  of  plotting 

D 

P.R.  vs  •  the  nominal  control  effort  (Jen)  and  picking  a  gain  that  maximizes 
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rxn  exists)  or  picking  a 


^P.R.  (in  the  case  when  no  specific  requirement  on  J 
gain  that  achieves  a  reasonable  trade  off  between  and  Jxn. 

Case  iii(b) ;  Jm  is  not  specified  (i.e.  no  concern  for  regulation  robustness). 

For  this  situation,  we  basically  design  for  robust  stability  as  discussed 
in  Section  III. 
j5_.  2_  Application  Example 

Consider  the  simple  scalar  system 

x  =  -x  +  u,  xQ  =  1  (5.6) 

min 

J  =  /  (x2  +  Pc  u2 )dt  (5.7) 

u  o 

Case  1 ;  Jm  is  given  and  e  is  given.  Let  k  =  3  so  that  Jm  =  3Jxn  and  let 

| A  a |  =  0.5  and  | A  b |  =  1.207  (5.8) 

Then  the  calculation  of  e,  pyp,  Jxn  an<i  Jm  is  summarized  in  Table  4. 

From  this  it  is  seen  that  e  is  greater  than  Myp  but  less  than  Uy  in  the 
range  of  pc  considered.  Thus  it  can  be  concluded  that  one  can  guarantee 
stability  for  the  given  perturbtions  but  cannot  guarantee  Jxp  <  kJxn  =  3Jxn. 
Either  the  range  of  perturbations  Aa,  Ab  have  to  be  smaller  to  be  able  to 
guarantee  Jxp  <  3Jxn  or,  for  the  given  Aa,  Ab,  one  has  to  sacrifice  more  in 
performance  by  allowing  higher  k  (i.e.  Jm).  For  example  if  k  is  increased  from 
3  to  37,  the  resulting  Pyp  and  Jm  are  as  shown  in  Table  5. 

For  this  situation,  it  is  seen  that  higher  the  control  effort  (lower  the 
pc)  better  is  the  performance  robustness  in  the  sense  that  higher  is  the 
perturbation  bound  for  robust  regulation  (Jxp  <  37Jxn). 

Case  ii(b)  :  Jm  is  given  and  e  is  not  given.  Let  Jm  =  JQ  (5.9a) 

oo 

where  JQ  =  open  loop  regulation  cost  =  j  x2  dt  (5.9b) 

•  o 

subject  x  =  -x  (5.10) 
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Let  JQ  =  kJxn,  where 

oo 

Jxn  ~  f  x2  dt  subject  to  (5.10). 
o 

Thus , 

Jo 

k(Pc)  = -  (5.11) 

Jxn  (Pc) 

indicating  that  k  and  Jxn  are  functions  of  the  control  weighting  Pc. 

The  table  of  the  perturbation  bound  Uyp  for  robust  regulation  is  given  by 
Table  6. 

Thus,  for  this  situation  also  it  is  seen  that  the  higher  the  control 
effort  (lower  the  pc)  better  is  the  performance  robustness  in  the  sense  that 
higher  is  the  perturbation  bound  for  robust  regulation  (JXp  <  JQ ) . 

Case  iii(a) ;  Jm  is  not  known  (k  is  not  given)  and  e  is  given  (and  regulation 
robustness  is  desired) . 

As  before  let  | A  a |  =  0.5,  | Ab |  =  1  .207.  This  time  Jm  is  the  design 

objective  variable  with  e  being  equated  to  hyp.  Table  7  summarizes  the 
different  quantities. 

Design  Observation  5.1:  As  anticipated,  lesser  the  e,  closer  is  k  to  1 . 

Design  Observation  5.2;  Comparing  the  entries  e  and  Jm  with  those  of  Table  4, 
it  is  seen  that  a  higher  Jm  is  to  be  tolerated  to  accommodate  the  given 
perturbation  Aa  and  Ab  of  (5.8). 

Design  Observation  5_.3_:  It  may  be  noted  that  in  this  case,  nominal  performance 
(Jxn)  and  performance  robustness  (^p.r.)  have  opposing  trends  indicating  the 

need  for  a  trade  off.  A  reasonable  choice  for  a  ’robust’  control  gain  could  be 

1  3 

the  gain  corresponding  to  Pc  =  -  or  -  or  1  •  However,  if  there  is  no  stringent 

2  4 

requirement  on  the  nominal  performance  then  the  control  gain  corresponding  to 
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higher  Pc  (such  as  2  above)  can  be  selected  since  that  gain  achieves  a  low 


value  for  Jm  for  the  given  perturbation. 


pc 

1/8 

1/4 

1/2 

3/4 

1 

5/4 

3/2 

2 

e 

2.914 

1 .992 

1  .38 

1.13 

1 

0.91  2 

0.85 

0.77 

pyp 

2 

1 .236 

0.732 

0.521 

0.414 

0.3416 

0.29 

0.224 

0.5 

0.6708 

0.865 

0.982 

1  .06 

1.118 

1.162 

1  .225 

3 

2.236 

1  .732 

1  .523 

1  .414 

1  .341 

1  .29 

1  .224 

TABLE  5 


pc 

1/8 

1/4 

1/2 

3/4 

1 

5/4 

3/2 

2 

e 

2.914 

1 .992 

1  .38 

1.13 

1 

0.912 

0.85 

0.77 

UYP 

2.918 

2.175 

1  .685 

1  .481 

1 .375 

1 .305 

1  .255 

1 .19 

6.18 

8.273 

10.68 

12.11 

13.08 

13.79 

14.33 

15.11 

Jxn 

0.167 

0.2236 

0.2886 

0.3273 

0.3535 

0.3726 

0.3875 

.4085 
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TABLE  6 


pc 

1/8 

1/4 

1/2 

3/4 

1 

5/4 

3/2 

2 

k 

3 

2.36 

1 .732 

1  .527 

1  .4149 

1.3416 

1  .29 

1  .224 

Myp 

2 

1  .236 

0.732 

0.527 

0.414 

0.3416 

0.29 

0.224 

0 

II 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

0.5 

3 

2.236 

1  .732 

1.5227 

1  .414 

1.3416 

1  .29 

1 .224 

Jxn 

0.167 

0.2236 

0.2886 

0.3273 

0.3535 

0.37265 

0.3875 

0.4085 

Pc 

1/8 

1/4 

1/2 

z 

2.914 

1  .992 

1  .38 

3 

2.236 

1  .732 

v  vy 

k= - 

Py-e 

34.84 

9.174 

4.926 

Jm=kJxn 

5.818 

2.05 

1  .42 

Jxn 

0.167 

0.2236 

0.289 

1 

®P.R.= 

0.172 

0.488 

0.704 

Jm 

— 

TABLE  7 


3/4 

1 

5/4 

3/2 

2 

1.13 

1 

0.912 

0.85 

0.77 

1  .525 

1  .414 

1.34 

1 .29 

1.224 

3.86 

3.42 

3.134 

2.932 

2.695 

1 .26 

1  .21 

1.16 

1.13 

1  .098 

0.327 

0.3535 

0.3726 

0.3875 

0.408 

0.793 

0.826 

0.862 

0.885 

0.91 
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VI.  EXTENSION  TO  LARGE  SPACE  STRUCTURE  MODELS 


Up  to  this  point,  the  analytical  development  to  analyze  and  synthesize 
robust  controllers  using  perturbation  bounds  was  presented.  Some  multiple 
input  multiple  output  system  examples  have  been  presented  to  illustrate  the 
concepts.  One  application  that  is  of  interest  to  the  present  research  under 
this  contract  is  that  of  the  control  of  Large  Space  Structures  (LSS).  It  is 
clear  that  the  proposed  development  can  be  easily  illustrated  using  an  LSS 
model  data.  However,  the  main  purpose  of  this  section  is  to  be  able  to 
exploit  the  special  structure  of  LSS  models  and  specialize  the  results  obtained 
so  far  to  the  case  of  LSS  models.  Towards  this  direction,  in  section  6.1  we 
obtain  explicit  expressions  for  the  perturbation  bounds  jj  ,  treating 
parameter  variation  as  perturbation,  in  terms  of  modal  data  and  draw  some 
conclusions  about  their  usefulness.  Then  in  section  6.2,  the  mode  truncation 
problem  is  treated  as  an  additive  perturbation  and  the  results  of  section  6.1 
are  applied  to  this  case  which  results  in  a  preliminary  algorithm  to  determine 
the  magnitude  of  the  control  gain  apriori  as  a  function  of  the  number  of  modes 
one  controls. 

6. 1  Uncertain  Modal  Data  as  Per turbat ion 

It  is  known  that  the  typical  differential  equation  describing  the  dynamics 
of  'modal  coordinates'  of  an  LSS  model,  assuming  no  control  inputs  (i.e.  open 

loop  case) ,  is  given  by 

••  . 

H  i  t  2  £  (i)  ^  +  (aj  =  0  (6.1) 

where  to  is  the  modal  frequency  of  ith  mode,  £  is  the  damping  ratio  assumed  to 
be  constant  for  all  modes  (0  <  £  <  1  )  and  ri^  is  the  ith  modal  coordinate.  The 
state  space  model  for  (6.1)  is  given  by 


-73- 


xi  -  Aii  xi 


(6.2) 


where 


xiT  =  [ni  nil 


(6.3) 


-i^±2  -2^0)1 


Clearly  the  system  is  stable  with  eigenvalues  -  £a)  ±  joj/1-  . 

We  are  now  interested  in  obtaining  perturbation  bounds  for  the  stability 
of  the  system  (6.2)  assuming  variations  in  the  matrix  A-^.  If  one  considers 
this  pertrubation  as  an  'unstructured*  perturbation,  which  amounts  to  ignoring 
the  structural  information  we  have  altogether,  the  elemental  bound  of  (2.1) 

can  be  obtained.  On  the  other  hand,  if  one  thinks  of  it  as  a  'structured' 
perturbation,  then  one  can  use  the  elemental  bound  PeY  °£  (2.5  )  using  Un 
matrix.  In  fact  one  can  do  even  better  by  not  only  considering  the 
perturbation  as  'structured'  perturbation  but  also  by  exploiting  the  additional 
information  that  the  variations  can  only  be  in  frequency  and  damping  ratio, 
thereby  using  the  matrix  Ue  having  the  structure 


0  0 


U 


(6.4) 


e 


These  three  bounds,  when  specialized  to  the  model  of  (6.3),  can  be  given 
explicitly  in  terms  of  the  modal  data  as  follows: 


The  Lyapunov  matrix  equation  is  given  by 


Ai±T  Pi  +  Pi  Aii  +  2I2  =  0 


(6.5) 


whose  solution  is  given  by 


-74- 


1  +  45  2 

Oli 

1  _ 

1/0^2 

+  —  ) 

2? 

1 

1  +  0)^2 

.2 

1 

(6.6) 


Note  that  in  this  case  P^m  —  (i.e.  all  the  elements  of  P-^  are  positive). 

The  bounds  UePi  (for  'unstructured'  perturbation),  ( 'structured '  but 

with  Un)  and  Uyi  ('structured'  with  Ue  of  (6.4))  are  given  by 


1  4ca)i3 

U  p.  =  -  [  -  - - - - ---] 


2  1  +  a^4  +  a)i2(2  +  4  ?2)  +  Si1/2 


1  8?oii3 

^eYi  =  "  f  ] 

2  Ui4  +  o)i2(2  +  4? 2 )  +  4?  aii  +  1  +  621/2 


(6.7a) 


(6.7b) 


where 


yi 


1  8?oii3 

=  _  [ - - ] 

2  a>  jL2  +  2?uii  +  1  +  631/2 


(6.7c) 


0!  =  1  +  a)i8  +  8  aii6  ?2  +  coi4  ( 1 65 4  -  2)  +  8  aii2?2 

02  =  (2  coi8  +  toi6  (4  +  16?2)  +  8 ?  aii5  +  u^4  (4  +  16?2  +  32^ 4 ) 

1 

+  Oil3  (16?  +  32? 3 )  +  a)i2  (4  +  1 6? 2)  +  8? aii  +  21  - 

4?  2aii6 


03  =  2uii4  +  aii2  (4  +  8? 2 )  +  2 


In  figure  6,  the  variation  of  the  bounds  with  respect  to  frequency  oi  (with 
constant  damping  ratio  ?  =  0.001)  is  presented.  From  this  figure,  it  is 
clearly  seen  that  both  bounds  ^eyi  and  Uy^  are  always  improved  bounds  over  Mg.pi 
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as  expected.  It  is  also  interesting  to  note  that  increases  without  bound 

as  the  frequency  increases  while  yeYi  reaches  a  maximum  (which  depends  on  £ )  at 
a  particular  frequency  and  then  decreases  as  a)  is  increased.  This  behavior  may 
have  some  useful  implications  in  synthesizing  robust  controllers.  In  fact  this 
aspect  is  seen  to  have  a  useful  extension  to  determine  the  trade  off  between 
the  number  of  modes  we  control  and  the  corresponding  control  gain  one  can  use 
to  achieve  a  given  amount  of  regulation  in  the  case  where  'mode  truncation*  is 
treated  as  an  additive  perturbation  and  only  control  spillover  occurs.  This 
topic  will  be  discussed  in  the  next  section. 

Finally,  the  variation  of  ajmax  (which  results  in  the  bound  Yi  with  Un) 
as  the  damping  ratio  £  is  varied  in  the  range  0  <  5  <  1  is  depicted  in  figure 
7.  From  this  figure  it  is  interesting  to  note  that  wmax  decreases  as  the 
damping  ratio  £  is  increased.  This  aspect  is  to  be  investigated  in  future 
research. 

6.2  Mode  Truncation  (Spi 11 over  Terms )  as  Perturbation 

One  design  philosophy  in  LSS  control  is  that  of  the  'Two-Model*  theory  of 
control  design  in  which  a  higher  order  'evaluation  model*  (which  serves  the 
role  of  the  physical  system  in  system  evaluation)  is  driven  by  a  lower  order 
controller  model.  In  this  framework,  the  modeling  error  is,  of  course,  'mode 
truncation* .  One  prominent  effect  of  mode  truncation  in  the  control  design 
process  is  the  instability  caused  due  to  the  interaction  with  residual  modes, 
labeled  'Spillover*  [56] .  There  has  been  a  considerable  amount  of  work  done  on 
compensation  for  or  elimination  of  'spillover* ,  notably  by  Balas  [56] ,  Sesak 
[57] ,  Lin  [58] ,  Longman  [59]  among  others.  In  this  section,  the  unmodeled 
dynamic  problem  is  treated  as  an  additive  perturbation  and  the  time  domain 
stability  robustness  analysis  reported  in  the  earlier  sections  is  extended  to 
the  case  of  spillover  instability.  Using  the  explicit  expressions  for  the 
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bounds  obtained  in  the  previous  section,  an  algorithm  is  proposed,  for  the 
case  of  control  spillover  alone,  by  which  one  can  determine,  a  priori  the 
magnitude  of  the  control  gain  in  relation  to  the  number  of  modes  one  wishes  to 


control . 

6.2.1  LSS  Models  and  Mode  Truncation  as  Perturbation 

Consider  the  standard  state  space  description  of  LSS  evaluation  model  with 
N  elastic  modes: 

X  =  Ax  +  Bu  +  Dw  x(0)  =  xo;  x+Rn=s2N#  u+Rm 


where 


y  =  Cx 
z  =  Mx  +  v 

yT  =  rv.T  v  T 


y+R* 


z+R* 


=  [x-|T,  X2T,...xnT];  = 


hi 


(6.8a) 


(6.8b) 


A  =  Block  diag.  A 


li 


0  1 


(6.8c) 


BT  =  [B1t,  B2T,...Bnt];  Bi  = 


J>iT_ 


( 6 . 8d  ) 


DT  =  [D1T,  D2T, . . .DNT] ;  Di  = 


-diT- 


(6.8e) 


C  =  [C-)C2.  .  .C^j]  and  M  =  [M-|  , M2 ,  .  .  . !%] 


w(t) 

w  o’ 

K 

(wT(t)  vT(x)]  = 

_v(t)_ 

_0  V 

(t-x) 

Let  the  above  system  be  evaluated  by  the  quadratic  performance  index 
1  im  1  t 

J  =  -  ?/  tyT(x )Qy(x )+uT(x )PcR0u(x ) ]dx 

t+a>  t  O 


(6.8f ) 

( 6 . 8g ) 


(6.9) 
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where  scalar  Pc  >  0  and  Q,R0  are  (kxk)  and  (mxm)  symmetric,  positive  definite 
matrices,  respectively  and  £  is  the  expectation  operator. 

Assuming  the  modes  n i  are  ordered  in  increasing  order  of  frequency,  let  us 
retain  the  first  Nr  modes  for  control  design  purposes.  Accordingly,  the 
reduced  order  control  design  model  of  dimension  nR(=2Nr)  <  n  is  given  by 

XR  —  ARxR  +  BRU  +  DRW  xr+R  k 

yR  =  CRXR  (6.10) 

Zr  =  MRXR  +  V 


where  the  above  control  design  model  is  obtained  by  direct  truncation  of  the 
full  order  model  given  by 


• 

XR 

Ar  0  ' 

XR 

br 

br 

X  = 

# 

= 

+ 

U  + 

_xT 

0  Aip 

Xiji 

Drp 

y  —  (CR  Ctp] 


XR 

XR 

z  =  [Mr  Mt] 

Xip 

Xip 

(6.11  ) 


Let  the  full  order  optimal  control  for  the  reduced  order  model  be  given 
by 


u  =  Gr  xr 


xr  =  ar  xr  +  +  Fr(z  -  Mrxr) 

A  A  /N 

=  ARXR  +  Frz  where  AR  =  Ar+BrGr  -  FrMr 


(6.12) 


where  GR  and  FR  are  the  standard  'controller'  and  'estimator'  gain  matrices 
respectively  [19]  such  that  the  closed  loop  system  matrix  for  the  control 
design  model  given  by 
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ar  ~ 


ar  brgr 

A 

_FrMr  Ar  _ 


(6.13) 


is  asymptotically  stable  under  the  usual  assumptions  of  controllability  and 
obs  e rvabi 1 i ty . 

The  closed  loop  system  for  the  evaluation  model  is  obtained  by  forcing  the 
evaluation  model  with  the  controller  of  the  control  design  model.  Thus,  we 
have 

* 

w 


— 

XR 

ar 

brgr 

0 

XR 

dr 

0 

• 

A 

A 

A 

= 

frmr 

ar 

FrMt 

XR 

+ 

0 

fr 

Xrp 

• 

0 

btgr 

A.ip 

XT 

Dip 

0 

— XT- 

_ 

y 

CR 

0 

C«p 

__u 

JD 

gr 

0 

[xR  XR  XT]T 


.  e .  ,  X  -  AbXb  +  Djj  w,  Yb  =  Gbxb 


(6.14a) 


(6.14b) 


where  the  stability  of  the  matrix  is  dictated  by  the  spillover  terms  BTGR  and 

F  rMiji  . 

We  now  write  Ab  =  A-|  +  Eb  where, 


(6.15) 


i.e.,  the  spillover  terms  are  treated  as  an  additive  perturbation.  Let  us 
denote  the  matrix  E  as  the  'spillover  matrix'.  Note  that  A^.,  ar  and  thus  A-\ 
are  asymptotically  stable  matrices  with  AT  having  a  block  diagonal  structure. 

Note  that  instability  of  the  closed  loop  system  may  occur  only  if  both  the 
'control  spillover  term  B^Gr'  as  well  as  the  'observation  spillover  term  FrMt>' 


ar 

brgr 

0 

0 

0 

0 

frmr 

ar 

0 

0 

0 

Fj^Mrp 

_0 

0 

at 

0 

btgr 

0  __ 
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are  present.  If  either  only  control  spillover  term  alone  or  observation 
spillover  term  alone  is  present,  then  clearly  the  closed  loop  system  is  stable 
but  significant  performance  degradation  may  occur.  At  the  present  stage  of 
research,  our  analysis  will  be  restricted  to  the  case  of  'control 
spillover'  alone.  Accordingly  for  this  case  the  closed  loop  system  matrix  is 
written  as 


A-j  —  Ai  +  Ei 


where 


(6.16) 


and 


ar  brgr  0 

0 

0  0 


Ar+BrGr  0 


T— 


(6.17a) 


0  0 
0  0 


0  BTGR 


0 

0 

0 


(6.17b) 


We  now  recall  the  expressions  and  of  (6.7)  &  (6.6)  respectively  to 

be  applied  in  our  present  analysis  on  mode  truncation.  From  (6.6) ,  (6.7)  and 

the  figure  6  the  following  design  observation  can  be  made. 


Design  Observation  6.1:  An  important  property  to  note  is  that  crmax (Pj_Ue ) s 
decreases  and  thus  increases  with  increase  in  tui* 

We  now  apply  the  above  'Perturbation  Bound  Analysis'  to  the  case  of  matrix 
A*|  of  (6.16).  From  (6.17)  it  can  be  seen  that  the  maximum  modulus  element  (in 
the  worst  case  situation)  of  the  'spillover  matrix*  will  be 

em  =  b  m  g  (6.18) 
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where 


b  lBij  Imax'  9  “  lGR. . Imax 

13 

(Note  that  bm  can  be  easily  determined  for  a  given  LSS  model  whereas  g  is  not 
known  before  design). 

Applying  (2.6)  to  the  matrix  Ai  of  (6.16),  we  conclude  that  Ai  (formed  by 
forcing  the  evaluation  model  with  a  reduced  order  controller)  is  stable  if 


em  < 


amax^pmue^s 


where  P1  satisfies  P1A1  +  A1TP1  +  2Ia  =  0  (a=n+nR)  and 


(6.19) 


U, 


eb2 


(because  of  the  structure  of  E  matrix  of  (6.17)) 


and 


U, 


eb2 


Ue1 1  •  ue12 
ue21  *  Ue22 


where  U( 


ei] 


0  0 
1  1 


(6.20) 


i.e.,  Ue^2  contains  NT  =  (N-Nr)  number  of  Ueij  submatrices  because  of  the 
special  structure  of  B  matrix  of  (6.8).  Note  that 


p1 


PR 

0 


0 


T_ 


where  PR  is  2nRx2nR  matrix  and 


PT  =  Block  diag  [...Pi...]  i  =  1,...,nt  (6.21) 

and  Pi  is  a  2x2  submatrix  having  the  same  structure  of  Pi  as  (6.6). 

Using  the  structural  information  of  Ue,  the  design  observation 
(6.1)  and  (6.7)  ,  it  may  be  shown  (after  some  matrix  manipulation)  that  we  can 
guarantee  avoidance  of  instability  due  to  spillover  if 
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I 


-  1 1  1 1  Pyi  Wyi 

g  °r  I  |  GR  |  |  <  <  — —  for  some  given  i=1,2,...N  (6,22) 

Nbm  N<pbm 

where  is  given  by  (6.8)  and  NT  is  the  number  modes  truncated  (=N-i). 

Even  though  it  is  known  that  control  spillover  alone  would  not  make  the 
closed  loop  system  (6.17)  unstable,  the  result  given  in  (6.22)  is  useful  from  a 
design  point  of  view  since  this  condition  specifies  apriori  a  bound  on  the 
control  gain  one  can  use  for  maintaining  stability  and  its  relationship  to  the 
number  of  modes  one  wishes  to  control. 

Design  Observation  6.2:  Note  that  the  bound  on  the  control  gain  for  (robust) 
stability  increases  with  increase  in  the  number  of  modes  retained  for  control, 
which  is  consistent  with  the  Spillover1  phenomenon. 

We  now  apply  the  above  result  to  the  LSS  example  considered  in  [60] • 
Example  6.1  :  The  example  considered  is  the  Solar  Telescope  (SOT)  model 

discussed  in  Yousuff ,  Skelton  [60] •  For  our  purposes  we  consider  the  model 
with  the  first  eight  elastic  modes  ordered  in  increasing  order  of  frequency  as 
the  evaluation  model.  These  modal  frequencies  are  given  by 


Mode  # 

1 

2 

3  1 

4 

5 

6 

7 

8 

“i 

r/s 

0.914 

3.43 

3.65 

i 

10.2 

14.8 

53.26 

149.37 

153.43 

The  matrix  B  (m=8)  given  in  [60]  is  not  reproduced  here  for  lack  of  space. 

The  damping  ratio  £  =  0.001  is  for  all  modes.  For  this  model,  the  quantity  bm 
of  (6.18)  is  given  by 

bm  =  0.03124 

The  tabulation  of  ^yj^  for  i  =  1,2,...N=8  (from  equation  (6.7c))  for  each 
o)j_  is  as  follows: 
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0)i 

yyi 

^yi/Nbrn 

0.914 

0.689x10-3 

0 .027x1 0“^ 

3.63 

0.559x10-2 

0.0223 

3.65 

0.563x10-2 

0.0225 

10.82 

0.178x10-1 

0.071 

14.85 

0.245x10-1 

0.098 

53.86 

0.892x10-1 

0.356 

149.37 

0.247 

0.988 

163.43 

0.254 

1  .016 

The  entries  in  the  last  column  suggest  that  the  more  modes  we  control,  the 
higher  the  tolerable  gain  for  stability. 

The  main  advantage  of  the  method  is  that  the  bound  on  the  control  gain 
depends  only  on  the  open  loop  data  such  as  the  modal  frequency,  mode  shape 
slopes  at  actuator  location,  the  damping  ratio  5,  the  number  of  control  inputs 
and  the  number  of  modes  included  in  the  evaluation  model. 

Of  course,  the  case  which  deals  with  combination  of  control  and 
observation  spillover  is  more  useful  and  interesting  and  more  research  is 
needed  to  extend  these  concepts  to  this  more  complicated  case.  It  is  to  be 
noted  that  even  though  ref.  [68]  treats  this  case,  the  main  result  stated  is 
valid  only  for  the  case  of  control  spillover  alone  and  this  restriction  is 
inadvertently  missing  in  the  statement  of  the  main  result. 
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Figure  6:  Variation  of  Bounds  with  frequency. 
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VII.  CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FUTURE  RESEARCH 


7 . 1  Work  in  Retrospec t 

The  main  theme  of  the  research  under  the  WPAFB  contract  has  been 
to  analyze  and  synthesize  controllers  for  robust  stability  and  performance 
for  linear  time  invariant  systems  subject  to  linear  time  varying  perturbations. 
First  the  aspect  of  stability  robustness  alone  was  considered.  The  main 
contribution  of  the  research  in  this  aspect  is  the  development  of  perturbation 
bounds  for  robust  stability  for  structured  uncertainty.  The  proposed  bound  is 
such  that  it  garners  the  structural  information  about  the  nominal  as  well  as 
perturbation  matrices  into  a  unified  expression.  Another  important 
contribution  is  to  reduce  the  conservatism  of  the  bounds.  A  solution  is 
proposed  to  improve  these  bounds  by  employing  a  scaling  transformation  and  the 
improvement  of  the  proposed  method  over  other  existing  methods  is  illustrated 
with  the  help  of  some  realistic  aircraft  control  examples.  Next  the 
application  of  these  bounds  in  the  analysis  of  stability  of  interval  matrices 
is  illustrated  and  the  extension  of  reduction  in  conservatism  of  Lyapunov 
stability  condition  to  ultimate  boundedness  control  of  linear  systems  is 
demonstrated  convincingly.  Finally  a  control  design  algorithm  is  presented  to 
synthesize  controllers  for  robust  stability  using  the  developed  bounds. 

Then  the  aspect  of  performance  (regulation)  robustness  is  addressed  and 
perturbation  bounds  are  again  presented  to  guarantee  prescribed  degree  of 
regulation  for  structured  uncertainty.  Another  important  contribution  is  to 
propose  a  design  algorithm  for  linear  regulators  for  robust  regulation  (which 
in  turn  guarantees  stability)  using  the  proposed  'Perturbation  Bound  Analysis'. 
The  most  significant  feature  of  the  proposed  design  methodology  is  that, 
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unlike  some  other  methods  in  which  the  gain  depends  explicitly  on  the 
perturbation  matrices,  the  proposed  gain  basically  is  a  nominal  gain  but,  of 
course,  implicitly  depends  on  the  perturbation  matrices  as  the  robustness  of 
the  gain  is  investigated  using  the  perturbation  bounds#  The  advantages  of  this 
concept  are  that  i)  one  need  not  have  any  restrictions  placed  on  the 
perturbation  matrices  AA,  AB  (in  fact  one  exploits  the  structure  of  AA,  AB  in 
the  design  procedure),  ii)  one  need  not  invoke  the  assumption  that  the  system 
be  controllable,  observable  throughout  the  range  of  perturbations,  and  iii)  the 
simplicity  of  the  control  law  is  maintained  in  the  sense  that  the  control  law 
is  still  designed  using  the  technique  used  for  nominal  design.  The  limitations 
of  the  method  could  be  that  for  given  perturbations  it  may  not  be  possible  to 
come  up  with  a  simple  control  law  that  accommodates  the  entire  given  range 
(however,  it  was  possible  to  do  so  for  a  specific  range  that  was  considered  in 
a  realistic  situation  as  demonstrated  by  the  example  in  section  2.3).  One  way 
to  circumvent  this  limitation  is  to  be  able  to  use  more  design  variables  and 
other  methods  to  improve  the  bounds  as  much  as  we  can. 

Next,  attention  was  focused  on  the  application  of  these  results  to  Large 
Space  Structure  models.  Utilizing  the  special  structure  of  these  models  and 
extending  the  results  on  the  perturbation  bounds  with  mode  truncation  as 
perturbation  it  was  possible  to  determine  the  control  gain  bounds  for 
maintaining  stability  in  the  case  of  control  spillover  alone  and  the 
relationship  of  these  control  gains  to  the  number  of  modes  one  controls.  This 
information  will  be  extremely  useful  in  analysis  asd  well  as  design  situations 
of  LSS  control  problems. 

The  publications  listed  as  Refs.  [61-70]  are  the  result  of  this  study. 

As  it  normally  occurs,  another  result  of  this  study  is  that  many 
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interesting  research  topics  surfaced  for  further  investigation •  These  are 
presented  in  the  next  section. 

7 . 2  Avenues  for  Further  Research  Which  Need  the  Continued  Support  of  the 

Air  Force 

1 )  The  foremost  area  of  research  would  be  to  further  reduce  the 
conservatism  of  the  perturbation  bounds. 

2)  It  is  important  to  develop  necessary  and  sufficient  conditions  for 

the  existence  of  a  linear  controller  such  that  (or  Bp.R.)  is  positive  for 

given  perturbations. 

3)  An  area  of  research  would  be  to  extend  the  development  of  explicit 
bounds  for  structured  perturbation  to  time-invariant  perturbations  and  examine 
the  reduction  in  conservatism  that  can  be  achieved. 

4)  Another  area  of  interest  is  to  compare  the  proposed  'Perturbation 
Bound  Analysis'  approach  to  design  with  other  relevant  methods  like  the 
Guaranteed  Cost  Control  of  Chang  and  P^ng  [13]  and  the  'multimodel  theory'  of 
Ackermann. 

5)  It  is  also  of  interest  to  probe  the  relationship  between  the 
perturbation  bound  and  the  corresponding  degree  of  stability  measured  by  a,  the 
real  part  of  the  dominant  eigenvalue. 

6)  Another  aspect  for  future  research  would  be  to  extend  the 
'Perturbation  Bound  Analysis'  for  actuator-sensor  location  problems. 

7)  An  area  of  extreme  interest  would  be  to  use  the  perturbation  bounds 
as  a  criterion  in  model/controller  reduction  and  develop  an  algorithm  for  same 
and  compare  it  with  other  relevant  schemes. 

8)  One  foremost  area  of  research  would  be  to  extend  the  proposed 
analysis  and  design  methodology  to  the  case  of  combined  modeling  errors  such  as 
parameter  variation,  mode  truncation  and  possibly  nonlinearities. 
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9)  There  is  need  for  probing  into  the  comparison  and  contrast  of 
frequency  domain  results  and  the  proposed  time  domain  results. 

10)  Of  course,  it  is  always  instructive  to  apply  the  developed 
methodology  to  practical  applications  such  as  aircraft,  spacecraft  control 
problems  and  robotics. 
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Appendix  A 


Proof  of  Theorem  2.1;  Consider  x  =  (A+E)  x(t)  (A1  ) 

where  lE^j  |max  =  e  (scalar)  and  A  =  eUn,  where  Un  is  an  nxn  matrix  with 
Un  =  1  for  all  i,  j=1,2,...,n. 

ij 

Let  V(x)  =  xTp  x  >  0  be  the  Liapunov  function  for  the  system  in  (A1  )  where 
P  is  the  symmetric  positive  definite  solution  of 


A?p  +  pa  =  -2ln 

Then,  V(x)  =  -xT  2Inx  +  xT(ETP  +  PE)x 
Now 

Let  e  < 


(A2 ) 
(A3) 


1 


Omax ( pmun ) ; 


*  amax(pmA)s  <  1 

+  CTmax^I>E^s  <  1 

-*■  |  A  ( PE )  g  |  max  <  1 

+  Ai[(PE)s  -  In]  <  0 

^-In  +  <PE>s^  negative  definite 

*  t-2ln  +  EtP+PE]  is  negative  definite 

V(x)  of  (A3)  is  <  0  for  all  x 

+  (A+E)  of  ( A1 )  is  stable 
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Appendix  B 


Proof  of  Theorem  2.3:  Let  A  be  stable  and  Ag  be  negative  definite • 
Let  crmax(E)  <  crmin(As) 


*  amax(Es)  <  ^min^Ag) 

*  |  A  i  ( Es )  |  <  ~Amax(As) 

**  ^max^s^  <  ”^max(Ag) 

+  Amax(As  +  Es )  <0  (from  Rayleigh* s  quotient  formula  for  symmetric 

matrices,  Noble,  1970). 

As  +  ES  is  negative  definite 

*  ReA^  t  (A  +  E) ]  <0 

■*  (A  +  E)  is  stable 
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